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INTRODUCTION 


The  objectives  of  this  project  were: 

•  to  ascertain  whether  the  course  of  spontaneous  or  phytohemagglutinin  (PHA)-induced 
blastic  transformation  of  human  lymphocytes  in  vitro  is  affected  by  exposure  to  con¬ 
tinuous  wave  (CW)  or  pulsed  (PW)  2450  MHz  microwaves; 

•  to  compare  effects  obtained  by  CW  and  PW  microwave  exposure  at  the  same  average 
specific  absorption  rate  (SAR); 

•  to  compare  effects  obtained  by  microwave  heating  and  conventional  heating  needed 
to  increase  the  temperature  of  the  sample  by  0.5,  1.0,  1.5  and  2.0°C. 

•  to  ascertain  whether  chromosomal  aberrations  are  induced  under  the  conditions  of 
these  experiments. 

In  order  to  provide  reliable  quantitative  data  on  exposure,  a  system  with  automated 
dosimetry  was  developed,  and  tested  for  biocompatibility  with  human  lymphocyte  cultures. 
A  method  for  quantitation  of  lymphocyte  transformation  based  on  image  analysis  was 
developed  and  tested.  Experiments  on  effects  of  CW  or  PW  2450  MHz  exposure  at  heating 
and  non-heating  levels  were  carried  out.  Effects  of  conventional  heating  were  examined 
and  compared  with  effects  of  microwave  heating.  For  convenience  and  clarity  the  work 
performed  is  described  in  the  following  sequence:  (1)  exposure  system  and  dosimetry,  and 
(2)  studies  on  lymphocyte  cultures,  and  (3)  conclusions. 

EXPOSURE  SYSTEM  AND  DOSIMETRY 


Description  of  the  Exposure  System 


The  experiments  planned  in  this  project  necessitated  the  design  and  assembly  of  an 
exposure  system  that  would  meet  several  engineering,  dosimetric  and  biological  require¬ 
ments.  The  system  should  provide: 


i  For 
Jt  l  ~ 


(1)  an  environment  compatible  with  lymphocyte  culture  in  vitro  in  respect  to  relative  a(J 
humidity  (90-100%),  CO2  content  in  air  (5%),  and  temperature  (37°C),  as  well  as  the  A_ 

possibility  of  heating  and  equilibrating  the  temperature  of  the  samples  up  to  39° C  __ 

01 

using  microwave  or  conventional  heating; 


on/ 

ilty  Codes 
AVd : 1  and/or 
h'pf'oiai 
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(2)  the  possibility  of  exposing  samples  (preferably  involving  the  exposure  of  multiple  sam¬ 
ples  under  identical  conditions)  to  CW,  or  PW  2450  MHz  with  a  high  ratio  of  peak 
to  average  power,  and  the  possibility  of  simultaneous  sham  exposure  of  control  sam¬ 
ples  under  strictly  comparable  conditions,  where  the  only  difference  between  exposed 
and  sham  exposed  samples  would  consist  in  the  presence  (or  absence)  of  microwave 
absorption  in  the  sample; 

(3)  continuous  measurements  and  registration  of  the  temperature  of  exposed  and  sham 
exposed  samples  during  exposure  to  provide  a  temperature  profile  over  time  history 
of  the  samples; 

(4)  continuous  monitoring  of  the  wave  form  and  repetition  rate  in  the  case  of  PW  expo¬ 
sures. 

The  first  and  second  requirements  were  met  by  placing  two  shorted  rectangular  S-band 
waveguides  300  mm  long  (Omega  Laboratories  Models  6101)  in  a  conventional  CO2  tissue 
culture  incubator  (Forma  Scientific)  thermostated  at  37°C.  A  styrofoam  block  placed  in 
the  waveguide  serves  as  support  for  the  sample  holder,  and  ensures  that  the  sample  is 
centrally  placed  in  the  waveguide  in  the  same  location  in  successive  experiments.  A  small 
hole  drilled  in  the  waveguide  allows  the  introduction  of  a  minimally-perturbing  temperature 
probe.  The  hole  is  drilled  near  the  guide  wall,  at  the  center  of  the  shorter  dimension,  this  is 
a  low  electric  field  point  and  the  hole  produces  minimal  disturbance.  One  of  the  waveguides 
serves  for  sham  exposures,  the  other  one  is  connected  through  a  matched  coaxial  cable  to 
the  microwave  power  source.  This  consists  of  several  elements.  The  output  from  a  CW 
Hewlett-Packard  Model  8616A  oscillator  feeds  a  Hughes  Model  1177H  travelling  wave  tube 
amplifier.  Pulsed  or  amplitude  modulated  waveforms  can  be  obtained  through  the  use  of 
a  Hewlett-Packard  Model  8403A  p-i-n  modulator,  for  high  peak  power  pulsed  exposures, 
the  signal  can  be  further  amplified  with  a  1  kW  amplifier  (MCL,  La  Grange,  IL).  The 
waveguide  receives  its  power  through  an  isolator  and  a  dual  directional  coupler  that  allows 
forward  and  reflected  powers  to  be  measured  with  Hewlett-Packard  Model  432A  power 
meters  and  waveforms  to  be  monitored  using  a  Phillips  PM  240  oscilloscope.  The  amount 
of  power  supplied  to  the  waveguide  is  controlled  by  the  out;  ut  from  the  oscillator  and  the 
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voltage  of  the  1  kW  amplifier,  and  monitored  by  measurements  of  the  forward  and  reflected 
powers,  and  by  the  readings  on  the  oscilloscope.  The  main  function  of  the  oscilloscope  is, 
however,  to  monitor  waveforms.  These  means  of  control  and  monitoring  are  only  auxiliary, 
as  the  conditions  of  exposure  are  characterized  by  dosimetry  based  on  measurements  of 
temperature  and  of  SAR  within  the  sample,  as  specified  by  the  3rd  requirement  the  system 
should  meet. 

The  dosimetric  part  of  the  system  consists  of  two  minimally-perturbing  temperature 
probes,  connected  to  a  Hewlett-Packard  Model  59306 A  relay  activator,  which  is  itself  con¬ 
nected  to  a  Keithley  Model  192  digital  voltmeter  (DVM).  Both  the  relay  activator  and  the 
DVM  are  under  the  control  of  a  Hewlett-Packard  Model  86  desk-top  computer  through  the 
HPIB  (IEEE-4SS)  bus.  Initially  a  Narda  Model  801  IB  non-perturbing  double  temperature 
probe  was  used,  later  two  Vitek  Model  101  probes  were  substituted  for  the  Narda  probe. 
The  Narda  probe  requires  cumbersome  calibration  against  a  NBS  standard  and  introduc¬ 
tion  of  correction  factors  to  relate  the  current  to  temperature,  a  linear  relationship  being 
maintained  over  limited  temperature  ranges.  Vitek  probes  are  much  easier  to  calibrate, 
are  more  stable,  and  the  current-temperature  relationship  is  linear  over  the  temperature 
range  of  interest. 

Under  computer  control  the  temperatures  of  exposed  and  sham  exposed  samples  are 
recorded  sequentially  and  stored  in  memory  (records  are  maintained  by  keeping  a  library 
of  micro-floppy  disks).  Since  the  temperature  of  the  sham  exposed  sample  is  essentially 
constant  in  view  of  the  stability  of  the  temperature  of  the  conventional  incubator,  it  can  be 
monitored  less  frequently  than  the  temperature  of  the  exposed  samples.  For  conventional 
thermal  exposures  the  temperature  can  be  elevated  by  increasing  the  temperature  of  the 
incubator,  maintaining  control  samples  in  smother  CO2  incubator. 

The  sample  holder  has  to  be  a  tissue  culture  vessel  made  of  materials  non- toxic  to 
cells  in  vitro  and  transparent  to  2450  MHz  microwaves.  The  geometry  of  the  container  has 
to  allow  for  uniform  energy  deposition  within  the  sample.  These  considerations  limited 
the  choice  of  possible  sample  holders.  Previous  experience  with  tissue  culture  exposures  in 
a  S-band  waveguide  demonstrated  that  the  use  of  T-flasks  or  Petri  dishes  leads  to  signifi- 
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cant  nonuniformities  in  energy  deposition.  Two  types  of  sample  holders  were  tested,  botk 
rectangular  tissue  culture  Lab-Tek  (Miles  Laboratories)  chamber  slides,  a  two  chamber 
model  No4802,  and  a  four  chamber  model  No4804.  Analysis  of  SAR  distribution  based 
on  measurements  of  temperature  (see  below)  in  slide  chambers  containing  tissue  culture 
medium  demonstrated  that  the  energy  distribution  within  and  between  chambers  of  the1 
two  chamber  model  is  nonuniform,  the  difference  between  chambers  being  about  50%.  In 
the  four  chamber  model  measurements  demonstrated  that  differences  in  energy  deposition 
remain  within  10%  of  an  average  value.  Therefore,  the  four  chamber  model  was  selected 
for  use  in  experiments.  Three  of  the  chambers  contain  tissue  culture  medium  and  cells 
the  fourth  contains  medium  only.  A  hole  is  drilled  in  the  cover  of  the  chamber  and  a  tern 
perature  probe  is  introduced.  In  this  way  each  experiment  can  be  carried  out  in  triplicate 
on  simultaneously  exposed  samples  with  simultaneous  continuous  dosimetric  control. 

Dosimetry  and  Temperature-Time  Analysis 

During  an  experiment  the  temperature  in  the  exposed  and  the  sham-exposed  samples- 
is  recorded  at  regular  intervals,  with  a  minimum  interval  of  less  than  1  s.  The  temperature 
is  recorded  before  the  beginning  of,  during,  and  after  the  exposure.  The  “on”  and  “off' 
times  of  exposure  are  recorded  on  the  computer  by  the  operator  using  a  soft-key  interrup : 
capability.  At  the  conclusion  of  a  run  the  T(t),  i.e.,  temperature  (T)  versus  time  (f  ■ 
behavior  is  analyzed  to  determine  SAR. 

If  the  SAR  is  5  (\Ykg~1 )  and  the  specific  heat  of  the  sample  is  C(Jkg-1K_1 ),  then  the 
rate  of  heating  during  microwave  exposure,  for  small  temperature  variations  over  which  C 
can  be  assumed  to  be  constant,  is 


If,  in  the  absence  of  deliberate  heating,  a  sample  is  above  or  below  its  equilibrium  temper¬ 
ature,  Teq,  with  its  surroundings,  then  for  small  temperature  differences  from  equilibrium, 
its  natural  rate  of  temperature  change,  from  Newton’s  law  of  cooling,  is  of  the  form 

(f)  =  «(r„-n  (2, 
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where  a  is  a  constant  that  depends  on  sample  geometry,  insulation,  etc.  In  general,  when 
microwave  power  is  applied  to  a  sample  that  was  not  in  thermal  equilibrium  with  its 
surroundings  at  the  start  of  exposure: 

§ <3> 

The  solution  to  eqn.  (3)  can  be  written  in  the  form 

r.T„=(^  +  T„_r„)(1_e-)  (4) 

where  T0  is  the  control  temperature  at  t  =  0.  Either  eqn.  (3)  or  eqn.  (4)  can  be  used  to 
determine  the  SAR.  From  eqn.  (3), 


where  The  subscript  5  indicates  the  application  of  microwave  power  corresponding  to 
SAR  5.  Therefore,  determination  of  the  rate  of  temperature  change  before  and  during 
application  of  microwave  power,  or  during  and  after  the  application  of  microwave  power, 
determines  the  SAR. 

If  the  sample  is  in  thermal  equilibrium  before  exposure  begins,  eqn.  (4)  reduces  to 
the  simpler  form 

T-T0  =  -e-°‘)  (6) 

Ql 

which  can  be  used  to  find  S  from  T(t).  In  general,  the  approach  using  eqn.  (5)  is  preferable, 
since  for  times  near  t  =  0,  the  exponential  cam  be  expanded  in  quadratic  form. 

In  practice,  this  analysis  is  performed  automatically.  At  the  end  of  a  run,  the  point 
of  inflection,  or  turning  point,  in  T(t)  is  found  numerically.  This  is  done  by  checking  the 
change  in  slope  of  the  temperature  as  a  function  of  time.  In  practice,  this  change  in  slope 
occurs  immediately  upon  the  application  or  removal  of  microwave  power,  within  the  time 
resolution  of  the  data  acquisition  system.  The  temperature  profile  to  the  left  and  right 
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of  this  point  is  fitted  by  least-squares  to  a  linear  or  quadratic  function  and  the  change  in 
slope  at  the  beginning  of  the  exposure  period  yields  the  SAR. 

In  practice  temperatures  in  the  exposed  sample  Eire  recorded  each  1,  2  or  3  seconds, 
the  readings  from  the  sham-exposed  sample  being  taken  at  10  times  longer  intervals.  In 
experiments  in  which  the  sample  is  microwave  heated  for  longer  periods  at  a  predetermined 
elevated  temperature,  measurements  are  made  less  frequently  once  equilibrium  is  reached, 
and  serve  only  as  control  of  the  stability  of  exposure  conditions. 

Measurements  made  in  air  or  water  in  the  incubator  outside  the  waveguides,  demon¬ 
strated  oscillations  of  ±0.2°C  resulting  from  the  on-off  cycle  of  heating  controlled  by  the 
thermostat.  Over  a  24  hour  period  such  oscillations  may  reach  ±0.5°C.  Measurements 
in  medium  in  chamber  slides  placed  within  the  waveguides  demonstrated  that  these  os¬ 
cillations  are  attenuated,  the  metal  structure  of  the  waveguide  acting  as  a  heat  sink  or 
source.  Over  a  24  hour  period  the  oscillations  do  not  exceed  0.05°C,  and  usually  remain 
within  0.02°C.  The  sensitivity  of  the  temperature  probes  is  nominally  0.01°C,  however 
compaxison  of  successive  readings  permits  extrapolation  to  0.005°C.  Based  on  data  in  the 
literature  the  biological  endpoints  examined  are  not  affected  by  temperature  excursions 
within  0.5°C. 

The  software  written  for  reading  temperatures  from  two  Narda  probes  and  for  “smooth¬ 
ing  and  plotting  of  experimental  data”,  j.e.,  for  SAR  computation  is  enclosed  as  attache- 
ment  1.  The  software  was  slightly  modified  to  accommodate  the  alternate  possibility  for 
taking  reading  from  Vitek  probes. 

The  dosimetric  system  described  above  can  be  used  in  experiments  on  the  effects  of 
heating  caused  not  only  by  microwaves  but  also  by  any  means  (eg.  conventional,  ultra¬ 
sound).  Because  of  its  wide  applicability  the  system  is  the  subject  of  two  brief  communi¬ 
cations:  one  at  the  IRPA-7  Congress  in  April,  1988  in  Sydney  Australia  (attachment  2), 
and  a  second  one  at  the  10th  Annual  Meeting  of  the  BEMS  in  June,  1988  in  Stamford, 
CT  (attachment  3).  Finally,  a  detailed  description  was  published  in  Health  Physics  56, 
303-307,  1989,  (attachment  4). 
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STUDIES  OF  LYMPHOCYTE  CULTURES 


Collection  of  Blood  Samples  and  Separation  of  Lymphocytes 

Approval  for  the  use  of  human  blood  was  obtained  November  19,  1986  from  the  Human 
Subjects  Research  Review  Board,  Office  of  the  Surgeon  General,  Department  of  the  Army. 
Volunteers  were  informed  about  the  nature  of  the  study,  potential  risks  and  the  potential 
benefit  derived  from  the  study,  and  signed  the  informed  consent  form,  as  approved  by  the 
HSRRR.  20ml  of  blood  was  withdrawn  by  venipuncture  under  aseptic  conditions  into  two 
heparinized  (143  USP  units)  sterile  evacuated  blood  collection  tubes  (Vacutaincr  No.  6480, 
Becton  and  Dickinson).  The  blood  was  transported  to  the  laboratory  at  room  temperature. 
Within  one  hour  from  withdrawal  mononuclear  cells  (lymphocytes  and  monocytes)  were 
separated  under  sterile  conditions  from  the  whole  blood  by  gradient  centrifugation  in 
Lymphocyte  Separation  Medium  (Boehringer,  Mannheim,  West  Germany)  according  to  the 
manufacturer's  instructions.  Following  separation  the  cell  suspension  was  washed  twice  in 
RPMI  1640  cell  culture  medium  (Gibco).  An  aliquot  was  withdrawn  for  the  determination 
of  the  cell  concentration,  cell  viability  and  cytological  (cytocentrifuge)  preparations,  the 
remainder  was  used  for  establishing  cultures. 

The  cells  were  counted  in  a  bright  line  hemacytometer  (Neubauer  type,  Spencer).  Vi¬ 
ability  was  tested  by  the  Trypan  blue  exclusion  test,  neutral  red  stain  for  the  vacuome,  and 
Janus  green  stain  for  mitochondria.  Cytocentrifuge  (Cytospin  2  cytocentrifuge,  Shandon) 
preparations  were  air  dried,  fixed  in  methanol,  stained  by  a  combined  Wright-Giemsa  stain 
and  cell  morphology  was  analyzed  under  a  microscope  (see  below). 

The  usual  yield  of  cells  ranged  from  lxlO7  to  3xl07.  The  suspension  consisted  of 
90-95%  lymphocytes,  the  remainder  of  cells  being  monocytes  and  occasional  (below  1%) 
granulocytes.  The  viability  of  cells  ranged  from  96.5  to  99.5%,  and  was  usually  about  98%. 

Cell  Cultures 

Four-chamber  culture  chamber  slides  were  filled  under  sterile  conditions  with  1  ml  of 
chromosome  medium  1A  (Gibco)  with  or  without  phytohemagglutinin  (PHA)  depending 
on  the  protocol  appropriate  for  the  particular  experiment.  106  cells  were  added  to  3  of  the 
chambers,  the  4th  contained  medium  only.  A  0.2  mm  hole  was  drilled  over  the  chamber 
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that  does  not  contain  cells  for  the  introduction  of  the  temperature  probe.  (For  experiments 
on  spontaneous  transformation  a  set  of  experimental  samples  consists  of  four  slides:  one 
with  PHA  (5  or  10  mg/m  1)  positive  control  incubated  in  a  conventional  incubator,  and 
three  without  PHA.  Out  of  these  one  is  incubated  in  a  conventional  incubator  at  37°C. 
one  is  microwave-exposed  and  one  is  sham  exposed). 

Cell  Harvest 

Following  incubation  the  cultures  without  PHA  were  mixed  within  the  culture  cham¬ 
ber  using  a  Pasteur  pipette,  and  transferred  to  microfuge  tubes.  The  cell  count  was 
established,  viability  was  tested  and  cytocentrifuge  preparation  were  made  (as  described 
above). 

In  the  case  of  cultures  incubated  with  PHA,  lOpl  of  colcemide  solution  (10  mg/ml, 
Gibco)  were  added  to  each  chamber  2  hours  before  harvest.  At  harvest  the  cultures  were 
transferred  to  microfuge  tubes,  spun  down  at  103  rpm  for  8  min  in  an  Eppendorf  micro  cen¬ 
trifuge,  and  the  supernatant  was  discarded.  1  ml  1%  sodium  cirtrate  solution  prewarmed 
to  37° C  was  added  with  vortexing.  The  cells  were  kept  in  this  solution  for  8  min,  spun 
down  under  the  same  conditions  as  above,  and  the  supernatant  was  discarded.  Fixative 
(cold  glacial  acetic  acid:  methanol,  1:3  vol/vol)  was  added  dropwise  with  vortexing.  After 
at  least  20  min  the  fixative  was  changed  once  or  twice.  A  droplet  of  the  cell  solution  in 
fixative  was  placed  on  a  microscope  slide  wet  with  20%  ethanol,  and  the  slides  were  flamed. 
Following  this  the  slides  were  air  dried,  and  stained  with  Giemsa  solution  in  phosphate 
buffer  pH  6.8  for  8  min  or  treated  with  trypsin  and  stained  with  Giemsa  according  to  the 
Seabright  procedure  (G-banding). 

Analysis  of  Microscope  Preparations 

Cytocentrifuge  slides  were  examined  under  a  100  x  oil  immersion  objective  and  the 
cell  and  nucleus  areas  were  determined  using  an  “Optomax”  image  analyzer.  The  original 
software  supplied  by  the  manufacturer  was  rewritten.  The  slides  was  coded  and  the  person 
who  made  the  measurements  was  not  aware  of  the  treatment  (exposed,  sham-exposed  or 
conventional  incubator).  The  code  number  of  the  slide,  the  area  of  the  cell,  the  area  of  the 
nucleus,  the  nucleus-cytoplasmic  ratio  and  a  symbol  for  morphological  features  of  the  cell 
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(S  for  a  small  lymphocyte,  1  for  an  intermediate  form  and  B  for  a  lymphoblastoid  cell) 
were  recorded,  and  stored  on  a  floppy  disk. 

The  process  of  lymphoblastoid  transformation  was  initially  evaluated  by  most  re¬ 
searchers  on  the  basis  of  cell  morphology.  Following  stimulation  the  lymphocyte  under¬ 
goes  gradual  enlargement,  accompanied  by  structural  alterations  of  the  nuclear  chromatin, 
and  becomes  an  intermediate  cell.  The  transition  of  the  intermediate  into  a  lymphoblas¬ 
toid  cell  consists  in  further  enlargement  of  the  cell,  accompanied  by  additional  changes  in 
chromatin  structure  and  the  appearance  of  nuclei.  Following  a  phase  of  DNA  synthesis 
(S-phase)  and  pre- mitotic  period  during  which  no  morphological  changes  occur  (g2  phase), 
the  lymphoblastoids  undergo  cell  division.  Evaluation  of  the  transformation  process  based 
on  qualitative  assessment  of  cell  morphology  is  subjective,  and  varies  with  different  ob¬ 
servers.  Because  of  this,  determination  of  tritiated  thymidine  incorporation  into  DNA  is  the 
most  widely  used  method  for  a  quantitative  assessment  of  lymphoblastoid  transformation. 
However,  this  method  does  not  fully  reflect  the  process  of  lymphoblastoid  transformation, 
being  only  a  measure  of  the  rate  of  DNA  synthesis.  Small  and  intermediate  lymphocytes, 
as  well  as  lymphoblastoids  in  other  stages  of  the  cell  cycle  rather  than  the  S-phase  are  not 
taken  into  account.  However,  using  this  method  it  is  not  possible  to  distinguish  between 
an  increase  or  a  decrease  in  tritiated  thymidine  incorporation  caused  by  a  change  in  the 
number  of  cells  synthesizing  DNA  or  by  a  change  in  the  rate  of  synthesis  alone.  As  the 
present  study  requires  incubation  of  cells  at  elevated  temperatures,  which  may  affect  not 
only  the  basal  metabolic  rate  but  also  the  DNA  synthesis  rate,  a  morphological  approach 
supplemented  by  objective  measurements  of  the  cell  and  nucleus  area  was  adopted. 

To  validate  this  approach,  three  different  observers  analyzed  cytocentrifuge  prepara¬ 
tions  from  cultures  incubated  with  PH  A  at  three  concentrations  of  5,  7.5,  and  10  mg/ml. 
The  last  one  was  indicated  by  the  manufacturer  as  the  optimal  concentration.  The  distri¬ 
bution  of  the  size  (area)  of  the  cells  demonstrated  a  dependence  on  PHA  concentration. 
The  area  of  the  nucleus  was  distributed  randomly  and  did  not  correlate  with  the  morpho¬ 
logical  classification  of  the  cell.  Consequently  the  nucleus-cytoplasmic  ratio  could  not  be 
used  for  the  evaluation  of  the  transformation  process. 
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Cells  of  size  between  80  and  120  pm2  were  classified  by  different  observers  as  small  or 
intermediate  lymphocytes,  and  cells  between  190  and  210  pm2  as  intermediate  or  blast' oid 
forms.  Therefore  based  on  morphology  and  measurements  of  the  cell  are?.,  cells  up  to 
100pm2  were  grouped  in  one  class  of  small  lymphocytes.  Cells  larger  than  200pm2  w/ere 
considered  as  blastoid  forms.  Intermediate  lymphocytes  were  additionally  subdivided  unto 
three  groups:  larger  than  100pm2  and  smaller  or  equal  to  120pm2,  larger  than  120 p.m2 
and  smaller  or  equal  to  150pm2,  and  larger  than  150pm2  and  smaller  or  equal  to  200p>:n2. 

The  distribution  of  the  cell  size  in  preparations  from  cultures  incubated  72h  w  ith 
PH  A  at  the  three  tested  concentration  is  shown  in  Table  1.  Significant  differences  were 
found  in  the  number  of  small  lymphocytes  (area  <  100/im2)  and  lymphoblastoids  (area  > 
200pm2).  The  decrease  in  the  number  of  small  lymphocytes  and  the  increase  in  the  numh.er 
of  blastoid  cells  with  increasing  PHA  concentration  validates  the  use  of  measurement:  -  of 
the  cell  area  for  the  assessment  of  lymphocyte  transformation. 

Table  1 

The  Distribution  of  cell  size  (area,  pm2)  following  72h  incubation 
with  PHA  at  Various  Concentrations  (mg/rnl) 


PHA  mg/ml 
Area  (a)  pm2 

5.0 

7.5 

10.0 

a  <  100 

42.5  ±  3.7 

26.8  ±  4.5 

17.4  ±  2.5 

100  <  a  <  120 

13.4  ±  2.0 

13.0  ±  2.8 

5.3  ±  1.1 

120  <  a  <  150 

15.6  ±  3.4 

12.1  ±  3.1 

10.0  ±  2.1 

150  <  a  <  200 

19.0  ±  2.6 

24.3  ±  2.8 

22.6  ±  3.0 

a  >  200 

10.6  ±  4.8 

25  ±  3.5 

46.2  ±  4.5 

Chromosomal  preparations  were  examined  using  a  100 x  oil  immersion  objective. 
Metaphase  olates  were  analyzed  for  numerical  structural  aberrations. 

Results  were  analyzed  using  the  “Statgraphics”  program  by  Statistical  Graphics  Cor¬ 
poration.  Analysis  of  variance  (ANOVA)  was  used  to  compare  the  numbers  of  cells  in 
particular  classes,  distributions  were  compared  using  multiple  regression  methods.  Results 
were  considered  significant  at  the  p  <  0.05  level. 
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Effects  Of  Exposure  On  Spontaneous  Transformation 

Blood  was  obtained  from  healthy  donors  by  venipuncture,  lymphocytes  were  separated 
as  described  above,  and  transferred  to  Chromosome  Medium  1A  without  PHA  (Gibco). 
The  cell  concentration  was  adjusted  to  106/ml.  One  milliliter  of  cell  suspension  was  in¬ 
troduced  into  each  of  the  3  chambers  of  the  chamber-slide,  the  4th  chamber  contained 
medium  only,  and  contained  the  sensor  for  the  dosimetric  system.  Thus  each  experiment 
was  carried  out  in  triplicate,  and  usually  repeated  3  times  (a  total  of  9  samples),  some 
experiments  were  carried  out  twice  (6  samples). 

Each  experiment  comprised  three  chamber-slides,  one  exposed,  the  second  one  sham 
exposed,  and  the  third  one  incubated  in  a  separate  incubator  set  at  37° C,  serving  as 
control. 

The  experiments  were  carried  out  in  three  series: 

•  exposed  to  C\Y  2450  MHz  at  a  non-heating  level,  i.e.  at  an  average  SAR  of  about  1 
W  /kg  (  ranging  from  0.S  to  1.3  W/kg),  and  at  heating  levels  resulting  in  an  increase  in 
temperature  of  the  sample  by  0.5°C  (SAR  1.8-2. 3  W/kg),  1°C  (SAR  3. 5-4. 5  W/kg),  1.5°C 
(SAR  G.8-S.3  W/kg),  and  2°C  (9.S-12.3  W/kg);  the  temperature  of  t lie  sample  varied 
over  the  incubation  period  by  ±0.2°C;  the  range  of  SARs  is  indicated  taking  into  account 
the  experiment  ally  observed  variations  and  the  ±10%  non-uniformity  of  energy  absorption 
over  the  three  chambers  of  the  chamber  slide; 

•  exposed  to  PW  2450  MHz  at  the  same  non-heating  and  heating  levels  of  temperature  and 
SARs  (for  these  PW  exposures  the  microwave  pulse  width  was  1/rs  and  the  pulse  repetition 
frequency  varied  from  100  to  1000  pps  as  the  exposure  changed  from  non-heating  to  the 
highest  heating  level); 

•  incubated  at  the  temperatures  obtained  by  microwave  heating  (37.5°C,  38°C,  38.5°C, 
and  39°C)  in  an  incubator  with  adjusted  thermostat,  samples  incubated  at  30°C  serving 
as  controls  and  comparison  to  the  non-heating  microwave  level. 

There  were  no  significant  differences  between  sham-exposed  and  control  cultures,  and 
therefore  the  results  from  these  were  pooled.  Table  2  presents  the  distribution  of  cell  sizes 
in  these  cultures. 
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Table  2 

The  distribution  of  cell  sizes  in  control  and 
sham-exposed  cultures  incubated  for  120h 


Area  (a) 

Number  of  cells 

pm2 

meamfcs.d. 

a  <100 

53.7  ±  2.9 

100  <  a  <120 

30.9  ±  1.8 

120<  a  <150 

10.1  ±  1.0 

150<  a  <200 

2.1  ±  0.8 

a  >200 

1.2  ±  0.9 

Table  3  presents  the  distribution  of  cell  sizes  in  conventionally  heated  cultures.  There 
are  no  significant  differences  between  cultures  incubated  at  30°C  and  37.5oC.  At  3S°C  the 
number  of  small  lymphocytes  (area<  100pm2decreases,  and  the  number  of  cells  with  areas 
between  120pm2and  200pm2increases,  thus  indicating  an  enhancement  of  the  transfoma- 
tion  process.  At  3S.5°C  the  number  of  cells  with  an  area  between  100  and  120 pm2 decreases, 
and  the  number  of  cells  with  area  bet  ween  150  and  200/nn2increases.  No  blastoid  cells  were 
found.  This  may  be  interpreted  as  enhancement  of  transformation  (blastoid  cells  divided 
and  the  daughter  cells  are  smaller)  or  an  inhibition.  The  latter  interpretation  is  less  likely 
because  of  the  results  obtained  at  38°C  ,  and  the  lack  of  cell  destruction  in  the  preparations. 


Table  3.  The  distribution  of  cell  sizes  in  control  cultures 
incubated  at  37.5,  38,  38.5,  and  39°C  , 
with  conventional  heating 


Number  of  cells 

mean±  s.d. 

Incubation 

temperature 
Area  (a)  (im2 

37°C 

37.5°C 

38°C 

38.5°C 

39°C 

a<100 

55.2±3.S 

54.8±2.2 

39.6±3.0 

53.0±4.5 

88.7±2.2 

100<  a  <120 

32.4±2.9 

29.4±3.6 

29.6±7.1 

14.7±2.3 

9.1±1.9 

120<  a  <150 

8.8±2.0 

11.2±1.7 

23.2±3.4 

16.2±2.2 

1.6±1.4 

150<  a  <200 

1.6±0.9 

1.3±1.0 

3.4±1.0 

15.8±2.0 

a  >200 

1.0±0.9 

2.4±1.7 

4.4±1.0 

0 

In  all  experiments  described  above  the  total  count  of  cells  did  not  change  significantly, 
being  of  the  order  of  106/ml  with  variations  within  the  error  of  the  method.  Cell  viability 
also  did  not  change,  remaining  at  the  initial  level  of  94  to  98%.  At  39°C  the  number  of 
small  lymphocytes  increased  significantly,  while  the  number  of  lymphoblastoid  cells  was 
smaller  than  in  cultures  incubated  at  37,  37.5,  and  38° C  ,  the  difference,  however,  was 
not  statistically  significant.  Cellular  debris  was  present  in  the  preparations.  The  total 
cell  count  decreased  from  6.7  x  105  to  9.1  x  105  cells/ml.  Viability  of  cells  decreased,  and 
ranged  from  84  to  93%.  These  results  may  be  interpreted  as  indicative  of  heat  cell  killing. 

Table  4  presents  the  results  obtained  following  CW  2450  MHz  microwave  heating. 
No  statistically  significant  differences  were  noted  between  this  and  the  preceding  series, 
with  the  possible  exception  of  cultures  heated  to  39°C  .  At  this  temperature  signs  of  cell 
destruction  were  noted.  The  total  cell  count  ranged  from  6.4  x  105  to  8.5  x  105  cells/ml. 
Viability  varied  from  75.5  to  81.5%.  In  view  of  this,  no  reliable  determination  of  the 
distribution  of  cell  sizes  could  be  made. 
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Table  4.  The  distribution  of  cell  sizes  in  cultures  exposed  to 
CW  2450  MHz  microwaves  at  non-heating  (37°  C  )  and  heating  levels 


Number  of  cells 

mean  ±  s.d 

Incubation 

temperature 
Area  (a)  pm2 

37°C 

37.5°C 

38°  C 

38.5°C 

39°C 

a<100 

53.3±2.1 

40.0±3.4 

3D.4±3.0 

43.6±4.3 

partial  cell 

100<  a  <120 

31.7±3.0 

36. Oil. 9 

30.0±4.9 

15.2±2.7 

destruction, 

120<  a  <150 

10.2±2.4 

12.2±2.0 

22.8±2.8 

17.6±1.7 

measurements 

150<  a  <200 

1.7±0.7 

4.8±1.6 

2.1±0.9 

22.6±2.8 

unreliable 

a  >200 

1.0±0.9 

2.4±1.7 

_ 

5.7±1.7 

1.0±0.9 

Table  5  presents  the  results  obtained  following  PW  2450  MHz  microwave  heating. 
Exposure  at  the  non-heating  level  results  in  an  increase  in  the  number  of  intermediate 
forms  with  an  area  greater  than  150pm2and  smaller  or  equal  to  200pm2,  the  number  of 
blastoid  forms  increases  too.  These  increases  may  be  considered  real  in  view  of  the  results 
obtained  at  37.5  and  38°C  .  At  these  temperatures  a  statistically  significant  increase  in  the 
numbers  of  lymphoblastoid  cells  (area<200pm2)  occurs.  At  38.5°C  cell  destruction  renders 
the  determination  of  cell  size  unreliable,  and  no  measurements  were  made.  In  view  of  this, 
planned  experiments  at  39°C  were  not  performed,  no  reliable  results  were  expected. 
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Table  5.  The  distribution  of  cell  sizes  in  cultures  exposed  to 
PW  2450  MHz  microwaves  at  non-heating  (37°C  )  and  heating  levels 


Number  of  cells 

mean  ±  s.d 

Incubation 

temperature 
Area  (a)  pm2 

37°C 

37.5°C 

38°C 

38.5°C 

39°C 

a<100 

59.7±3.9 

44.6±3.8 

46.8±4.7 

partial 

100<  a  <120 

15.4±1.8 

20.6±2.4 

16.3±2.0 

cell 

120<  a  <150 

13.7±1.4 

1G.6±1.5 

16.3±2.0 

destruction 

not 

150<  a  <200 

4.9±1.S 

3.9±1.0 

10.6±2.0 

measurements 

performed 

a  >200 

6.3±1.7 

14.4±1.2 

15.9±3.7 

unreliable 

Effects  of  Exposure  on  PHA-induced  Transformation 

Blood  was  obtained  from  healthy  donors  by  venipuncture,  lymphocytes  were  sepa¬ 
rated  as  described  above,  and  transferred  to  Chromosome  Medium  1A  containing  5  mg/ml 
PHA  (Gibco).  The  cell  concentration  was  adjusted  to  106  cells/ml.  One  milliliter  of  cell 
suspension  was  introduced  into  each  of  3  chambers  of  the  chamber  slide,  the  4th  cham¬ 
ber  contained  medium  only,  and  housed  the  sensor  of  the  dosimetric  system.  Thus  each 
experiment  was  performed  in  triplicate,  and  usually  was  repeated  three  times:  in  some 
instances  twice.  Each  experiment  comprised  three  chamber-slides,  one  exposed,  a  second 
one  sham-exposed,  and  a  third  one  incubated  at  30°C  in  a  separate  incubator,  serving  as 
control.  The  conditions  chosen  for  these  experiments  were  based  on  previously  obtained 
results.  No  experiments  were  carried  out  under  conditions  that  did  not  demonstrate  signif¬ 
icant  differences  between  experimental  and  control  samples,  or  resulted  in  cell  destruction 
rendering  measurements  unreliable.  The  experiments  were  carried  out  in  three  series: 

•  exposed  to  CW  2450  MHz  microwaves  over  the  whole  72h  incubation  period  at  a  non¬ 
heating  level,  i.e.  at  an  average  SAR  of  about  lVV/kg  (ranging  from  0.7  to  1.3  W/kg) 
and  at  a  level  resulting  in  an  increase  in  the  temperature  of  the  sample  by  1°C  (38°C  ) 
i.e.  at  an  average  SAR  of  about  4W/kg;  the  temperature  of  the  sample  varied  over  the 
course  of  the  experiment  by  ±  0.2°C  ;  the  indicated  range  of  SAR  takes  into  account  the 
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experimentally  observed  variations  in  different  chambers  of  the  chamber-slide; 

•  exposed  over  the  whole  72h  of  the  incubation  period  to  PW  2450  MHz  microwaves  at 
the  same  levels  of  average  SAR  and  temperature  as  in  the  preceding  series;  an  additional 
set  of  experiments  consisted  in  exposure  to  PW  2450  MHz  at  an  average  SAR  resulting  in 
an  increase  in  the  temperature  of  the  sample  by  1°C  over  the  first,  second  and  third  24h 
of  the  72h  incubation  period,  following  which  the  cultures  were  harvested. 

•  cultures  incubated  for  72h  in  an  incubator  adjusted  to  38°C  and  compared  to  cultures 
incubated  at  37°C  . 

The  results  of  these  experiments  are  summarized  in  Table  G. 


Table  G.  The  distribution  of  cell  sizes  following  721i 
exposure  to  PW  2450  MHz  microwaves  at  a  non-heating  (37°C  ) 
and  heating  level  (3S°C  )  compared  to  sham-exposed  samples. 


Number  of  cells 

mean  ±  s.d 

Incubation 

Shams  (37°C  ) 

37°C 

3S°C 

temperature 

Area  (a)  f, ini 2 

a<100 

43.2±7.4 

14.3±5.1 

11.7±5.3 

100<  a  <120 

14.3±4.6 

16.7±3.0 

5.9T2.9 

120<  a  <150 

14.9±5.3 

20.0±G.2 

7.8±2.7 

150<  a  <200 

11.7±5.1 

19.9±4.6 

14.5T4.1 

a  >200 

15.9±6.4 

29.1±6.6 

57.7±10.4 

Control  and  sham-exposed  cultures  did  not  differ  significantly  from  each  other.  Heat¬ 
ing  to  38°C  by  conventional  means  and  by  CW  microwave  exposure  enhanced  transfor¬ 
mation  to  a  limited  extent.  Exposure  to  CW  microwaves  at  a  non-heating  level  did  not 
affect  the  transformation  process.  Exposure  to  PW  2450  MHz  during  the  whole  incubation 
period  of  72h  enhanced  transformation  at  both  the  heating  and  non-heating  levels.  Expo¬ 
sure  to  PW  2450  MHz  at  the  heating  (38°C  )  level  during  the  first  24h  of  the  incubation 
period  did  not  affect  transformation,  no  significant  differences  between  exposed  and  sham- 
exposed  samples  were  demonstrated.  Exposure  during  the  second  24h  of  the  incubation 
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period  resulted  in  an  inhibition  of  transformation  arid  a  complete  absence  of  cell  divisions. 
Exposure  during  the  last  24h  enhanced  transformation. 

No  chromosomal  aberrations  were  demonstrat'  1  under  any  exposure  conditions. 

CONCLUSIONS 

Spontaneous  transformation  of  human  lymphocytes  in  vitro  is  enhanced  by  conven¬ 
tional  heating  and  CW  2450  MHz  microwave  heating.  No  significant  differences  wrere 
demonstrated  between  preparations  obtained  from  cultures  heated  to  the  same  tempera¬ 
ture  by  conventional  means  and  by  CW  microwaves.  PW  microwave  exposure  enhances 
transformation  at  a  non-heating  level.  Further  enhancement  of  transformation  was  seen 
in  preparations  heated  by  PW  exposure  to  37.5  and  3S°C  ,  the  number  of  transforming 
cells  being  significantly  higher  than  in  cultures  heated  by  CW  exposure.  Following  PW 
exposure  at  the  heating  level  resulting  in  an  increase  in  the  termperature  of  the  sample 
by  1.5°C  (to  38.5°C  )  or  2°C  (to  39°C  )  numerous  cells  undergo  destruction  making  the 
interpretation  of  the  preparation  unreliable. 

Effects  of  CW  exposure  seem  to  be  related  to  heating,  i.e.  depend  upon  the  increase 
in  temperature.  Effects  obtained  following  PW  exposure  differ  from  those  seen  following 
CW  exposure,  PW  exposure  enhances  transformation  at  non-heating  levels,  while  CW 
exposure  does  not.  Thus  effects  of  CW  and  PW  2450  MHz  microwaves  are  different  and 
are  not  related  to  the  average  SAR. 

PHA-induced  transformation  is  affected  by  exposure  in  a  similar  manor,  and  leads  to 
the  same  conclusions.  Additionally  this  part  of  the  study  demonstrated  the  dependence  of 
the  effects  on  the  stage  of  the  transformation  process.  Exposure  during  initial  stages  (first 
24h  of  the  72h  incubation  period)  is  without  effect.  Exposure  during  the  second  24h  of  the 
72h  incubation  period  blocks  cell  division  and  decreases  transformation,  while  exposure 
during  the  last  24h  enhances  transformation. 

Results  of  this  study  were  presented  at  the  10th  Annual  Meeting  of  the  Bioelec¬ 
tromagnetics  Society  (2  presentations)  in  Stamford,  Connecticut  in  19S8  and  at  the  7th 
International  Congress  of  the  International  Radiation  Protection  Association  (IRPA)  in 
Sydney,  Australia  (2  presentations).  Copies  of  abstracts  are  included  in  this  report.  One 
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paper  describing  the  exposure  system  and  dosimetry  protocols  has  already  been  published 
in  Health  Physics  (copy  included).  Two  further  manuscripts  are  in  preparation  and  will  be 
submitted  later.  Also  included  as  an  attachment  to  this  report  is  the  computer  software 
used  for  experimental  monitoring,  temperature  history  data  acquisition,  and  dosimetry. 
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Attachment  #1 


PROGRAM  FOR  COLLECTING  DATA 


50 

50 

70 

BO 


20  DISP^ARpYOU  USING  VITEK  OR  NARDA  PROBES, TYPE  V  OR  N" 

30  INPUT  PRBTYPS  ^  ^  .  .  . 

i4i)  PR0BC1-9.6184  9  PR0BC2-41 .459  !  Calibration  factors  for  Naraa  probe  1 
PRGBC3-9.9121  5  PR0BC4-42.386  !  Calibration  factors  for  Narda  proDe  2 
IF  PRBTYP$-"N"  THEN  GOTO  90  !  Probe  type  Narda  -  don't  reset  calib.  factors 
PRDBCl-100  §  PR0BC2-0  !  Calibration  factors  for  Vitek  probe* 

PR0BC3-100  #  PRDBC4-0  !  Calibration  factors  for  Vitek  probe 
90  ON  KEY*  1. “POWER  ON”  G0SUB  1360 
100  ON  KEY*  4, -NEW  OBS  INT"  GOTO  1100 
110  ON  KEY*  7. “POWER  OFF*'  G0SU&  1390 

120  HASS  STORAGE  IS  ":D500"  „  nnrrT/C04 

130  DIH  ATEMP< 1 000) , ATIHE( 1000) .BTEMP( 1000) ,BTIME< 1000) , PONT (50) ,POFFT(50 ) 

T40  DIH  ANUH(20) ,BNUM(20) , I PON (20) ,IPOFF<20) 

160  !  BEEP  A0, 100  9  WAIT  50  9  BEEP  A0, 100  9  WAIT  50  9  BEEP  A0, ICO  9  WAIT  50  *  BE 
IP  80,150 

170  DISP  "ENTER  DAY .MONTH , YEAR" 

180  INPUT  EXPDATE5 
ISO  EXPDATE-0 

200  DISP  -ENTER  SAHPLE  IDENTIFIER" 

210  INPUT  SAHPLES 

220  EXPDATAS-EXPDATESASAMPLES 

230  DISP  "ENTER  T1HE  OF  DAY  -  2A  HOUR  CLOCK"  i  INPUT  TODS 

2A0  ON  ERROR  GOTO  300 

250  CREATE  EXPDATAS6" :D501" , 1  .32000 

260  ASSIGN*  1  TO  EXPDATA5&" :D50 1 " 

270  !  PRINT*  1  ;  EXPDATES 
280  !  PRINT*  1  ;  SAMPLES 
290  GOTO  390 

300  IF  ERRN  -128  THEN  DISP  "DISC  IS  FULL  -  INSERT  A  NEW  DISC  AND  START  AGAIN" 
310  IF  ERRN  -128  THEN  GOTO  3A0 
320  IF  ERRN  -S3  THEN  GOTO  350 

330  DISP  "AN  ERROR  HAS  OCCURRED  -  PRESS  CONTINUE  TO  START  AGAIN" 

3A0  PAUSE  *  GOTO  150 

350  DISP  "THIS  DATA  FILE  NAME  ALREADY  EXISTS  -  ENTER  ANOTHER  SAMPLE  NUMBER” 

360  INPUT  SAMPLES 

370  EXPDATAS-EXPDATES&SAMPLES 

380  GOTO  250 

390  KK-0  !  Counter  for  new  observation  times  for  probes  A  and  B. 

400  BIN-0  !  Holds  previous  number  of  observations  from  probe  A. 

410  NDPTS-Q  !  Counter  for  number  of  data  points  recorded. 

420  CLEAR 

430  DISP  "This  program  will  store  1000  temps,  and  times  for  both  proDes. 

440  DISP  -ENTER  TDTAL  OBSERVATION  TIME  IN  SECONDS." 

450  INPUT  OBSTIME 

460  DISP  "Probe  A  should  be  in  the  exposed/experimental  culture.  * 


470  DISP 

480  DISP  "ENTER  OBSERVATION  INTERVAL  FOR  PROBE  A  IN  SECONDS." 

490  INPUT  NA 

500  IF  NA>0BSTIHE  THEN  520 
510  GOTO  540 

520  DISP  "You  selected  an  obs.  interval  >  the  total  obs.  time;  try  again. 
530  GOTO  430 

540  AJOBS-INT  (OBSTIHE/NA) 

550  NA-NA-1000 
560  H0LD-AJ0BS+8IN 
570  IF  H0lD>1000  THEN  590 
580  GOTO  610 

590  DISP  "You  have  requested  too  many  readings  of  probe  A;  try  again." 

600  GOTO  420 


610  DISP  "The  observation  interval  for  probe  B  should  be  an  integer  multiple" 
620  DISP  "of,  and  ideally  <  or  •  10  times  that  for  probe  A." 

630  DISP 

640  DISP  "ENTER  OBSERVATION  INTERVAL  FOR  PROBE  B  IN  SECONDS." 

650  INPUT  NB 

660  BJOBS-INT  COBSTIME/NB) 

670  NB-NB-1000 

630  IF  FP  (NB/NA>#0  THEN  700 
690  GOTO  720 

700  DISP  "The  obs.  time  for  B  is  not  a  multiple  of  that  for  A;  try  again." 

710  GOTO  640 

720  REMOTE  502  f  DVM  to  remote. 

730  CLEAR  6  KEY  LABEL 

740  REMOTE  507  !  Relay  to  remote. 

750  LOCAL  LOCKOUT  5  !  Local  lockout  for  DVM  and  relay. 

760  OUTPUT  502  ;"F0R2T0X"  !  DVM  to  DCV.2V  range,  and  cont.  on  talk. 

770  OUTPUT  507  "A12”  !  Connect  all  A1&A2  terminals  to  C16C2  terminals. 

780  ENTER  502  ;  AS  !  Trial  read  of  DVM. 

790  IJ1-0  1  Set  counter  for  number  of  times  power  turned  on. 

800  IJ2-0  !  Set  counter  for  number  of  times  power  turned  off. 

810  J-1  !  Set  counter  for  reading  probe  B. 

820  NT-NB/NA  ?  Set  interval  between  reading  probe  B. 

830  IF  KK> 0  THEN  860 

840  SETTIME  0,0  !  Zero  clock  for  easy  interpretation  only. 

850  !  Loop,  read  probe  A,  read  time,  put  values  into  arrays. 

860  NNA-1  «  NNB= 1 

870  FOR  1-1  TO  AJOBS  £  ENTER  502  ;  ASS  ATIME(I)-TIME  +<DATE  -EXPDATE>*86400  £  A 
£MP<I)-VAL  ( AS [ 5 ] ) 

880  ATEMP( I )-PR0BC1 »ATEMP< I )+PR0BC2  !  Evaluate  temp,  from  voltage  reading. 

890  NNA-NNA+1 

900  DISP  "A TEMP-" , ATEMP< I> ,"TIME-",ATIME(I> 

910  IF  I-J»NT  THEN  930  !  Check  to  see  if  probe  B  is  to  be  Tead. 

920  GOTO  1010 

930  OUTPUT  507  "B1”  !  Switch  B1  to  Cl  prior  to  reading  probe  B. 

940  WAIT  200  !  Wait  200ms  for  DVM  to  read  probe  B. 

950  !  Read  probe  B,  read  time,  put  values  into  arrays.  Switch  A1  to  Cl. 

960  ENTER  502  ;  B$6  BTIME< J)-TIME  ♦(DATE  -EXPDATE >*»86400  6  BTEMP(J>-VAL  (BSI53) 
£  OUTPUT  507  "A  1 " 

970  8TEMP( J)-PR0BC3»BTEMP< J)+PR0BC4  !  Evaluate  temp,  from  voltage  reading. 

980  NNB-NNB'-H 

?5SnDISP  "BTEhp“M.BTEMP<J>,"TlHE-",BTIME<J>  $  DISP  8  DISP  £  DISP  £  KEY  LABEL 
1000  J-J+1  f  Increment  counter  on  read  probe  B. 

1010  WAIT  NA  !  Wait  time  between  reading  probe  A. 

1020  NEXT  I  ♦  Increment  loop. 

1030  LOCAL  5 

Incn  "B12"  !  Simply  turn  off  lights  on  relay. 

1060  DISP  #^DISPn£  ^a*p  °n  £:*iSC•  *ou  be  Prompted  for  more  inputs  shortly." 

1070  WAIT  2000 
1080  GOSUB  DISCPRINT 
1090  GOTO  1400 
1100  CLEAR 


nrcD  ”5°  pr*ss  any  keys  until  the  computer  responds  with  a  beep  and" 
displays  Data  storage  is  complete'.  You  will  then  be  prompted  to"" 


1110  DISP  "Do  not  press  any 
1120  DISP  "displays  'Data  s 
1  30  DISP  "key  in^new  obse? 
1140  WAIT  5000 
1150  GOSUB  DISCPRINT 
1160  GOTO  440 

ii&a  ^J^CPRINT ‘  f  Write  all 
1 180  PRINT#  1  ;  NNA- 1 
1190  FOR  1-1  TO  NNA-1 


new  observation  data.' 


data 


disc 


this  observation  interval. 


1210 

1220 

1230 

1240 

1250 

1260 

1270 

1280 

1290 

1300 

1310 

1320 

1330 

1340 

1350 

1360 

1370 

1380 

1390 

1400 

1410 

1420 

1430 

1440 

1450 

1460 


NEXT  I 

PRINT*  1  ;  NNB-1 
FOR  1-1  TO  NNB-1 
PRINT*  1  ;  BTEMP( I )  ,BTIM£ <  I  > 

NEXT  I 

PRINT*  1  ;  IJ1  !  The  number  of  times  power  turned  on, 
IF  IJ1-0  THEN  1290 
FOR  1-1  TO  IJ1  9  PRINT*  l  ; 

PRINT*  1  :  IJ2  !  The  number 
IF  IJ2-0  THEN  1320 
FOR  1-1  TO  IJ2  $  PRINT*  1  ; 


PONT < I >  9  NEXT 
of  times  power 


turned  of f . 


POFFT ( 1 > 


TODS  !  Writes  time  of  day 


9  NEXT  I 
to  disc 


complete."  §  RETURN 


PRINT*  1 
KK-KK+ 1 
BIN-AJDBS 

NDPTS«NDPTS+2*(NNA-*-NN6)+IJ!  +  IJ2*lO 
BEEP  9  CLEAR  6  DISP  "Data  storage  is 
GOTO  1400  _  „ 

BEEP  9  IJ1-IJ1+1  9  DISP  "Power  on"  9  P0NT( I J 1 ) -TIME  9  RETURN 
BEEP  6  IJ2-IJ2+1  6  DISP  "Power  off"  6  POFFT ( I J21-TIME  9  RETURN 
DISP  "Do  you  wish  to  record  more  data7  ENTER  (Y)ES  or  (N)D." 

INPUT  REPLYS 
IF  REPLYS«"Y"  THEN  420 
ASSIGN*  1  TO  * 

DISP  "The  final  disc  management  is  now  being  completed" 

ASSIGN*  1  TO  EXPDATASA" :  D50 1 "  „  „  _  ,nrr_ 

IA1-1  9  IATOT-O  9  IB1-1  9  IBTOT-O  9  ION! - 1  9  IONT-O  e  I0FF1-1  «•  iOFFT- 


1470  FOR  J1-1  TO  KK 

1480  READ*  1  ;  ANUM(J1>9  IAT0T-IAT0T+ANUM( Jl ) 

1490  FOR  I-IA1  TO  IATOT  9  READ*  1  ;  ATEMP(I) ,ATIME<I>9  NEXT  I 
1500  READ*  1  ;  BNUM(J1>9  IBT0T-IBT0T+BNUM< J1 > 

1510  FOR  I-IB1  TO  IBTOT  9  READ*  1  ;  BTEHPU)  ,BTIME<I >9  NEXT  I 

toZU  RtAu*  i  ;  ii-uNiji>9  iGn  i -iUN  i  firUNt  jl )  9  IF  IPON(J1>-0  THEN  1540 

1530  FOR  I-I0N1  TO  ICJNI  <*  KtHu*  i  *.  ruin  (i  it  Nca.  I 

1540  READ*  1  ;  IPOFF(J1)C-  IGFFT-IuFFT  +  IPOFFt  Jl )  9  IF  IPOFF(J1>-0  THEN  1560 

1550  FOR  I-I0FF1  TO  IOFFT  9  READ*  1  ;  P0FFT<I>9  NEXT  I 

1560  I A 1 -IA1 +IAT0T  9  IB  1  - IB 1 ♦ IBTOT  9  I0N1 -I0N1 +I0NT  9  I0FF1 -I0FF1 +I0FFT 

1570  NEXT  Jl 

1580  READ*  1  ;  TODS 

1590  ASSIGN*  1  Tu  - 

1600  PURGE  EXPDATASA" :D50 1 " 

1610  PACK  ":D501" 

1620  NDPTS-8-NDPTS 

1630  CREATE  EXPDATA44":D501",1 .NDPTS 
1640  ASSIGN*  l  TO  EXP0ATASA" : D50 I" 

1650  PRINT*  1  ;  IATOT  9  DISP  "TOTAL  PROBE  A  READINGS  -  ", IATOT 

1660  FOR  1-1  TO  IATOT  9  PRINT*  1  ;  ATEMP(I) ,ATIME<I )  9  NEXT  I 

1670  PRINT*  1  ;  IBTOT  9  DISP  "TOTAL  PROBE  B  READINGS  -  ", IBTOT 

1680  FOR  1-1  TO  IBTOT  9  PRINT*  1  ;  BTEMP<I> ,BTIME(I)  9  NEXT  I 

1690  PRINT*  1  ;  IONT  9  DISP  "MW  POWER  TURNED  ON  ",I0NT,"  TIMES"  9  FOR  1-1  TO  101 

T  9  PRINT*  1  ;  P0NT(I>  9  DISP  P0NT<I>  9  NEXT  I 

1700  PRINT*  1  ;  IOFFT  9  DISP  "MW  POWER  TURNED  OFF  ", IOFFT,"  TIMES"  9  FOR  1-1  TO 
IOFFT  9  FRINT*  1  ;  POFFT(I)  9  DISP  P0FFT(I>  9  NEXT  I 
1710  PRIM*  1  •  TODS  9  DISP  "Time  of  day", TODS 
1720  ASS  I  IN*  1  TO  * 

1730  DISP  "Sample  number  is  ".SAMPLES 
1740  DISP  "Date  of  exposure  is  ".EXPDATES 

1750  DISP  "Finished  the  disc  management.  GOODBYE,  HAVE  A  NICE  DAY!" 

’760  END 


PROGRAM  FOR  DATA  ANALYSIS ,  CURVE  SMOOTHING  AND  GRAPHICS  version  A 


;o  ! . . . . * . ** 

20  !  * 

30  !  -  PROGRAM  TO  READ  TEMPERATURE  READINGS  TAKEN  FROM  * 

AO  »  *  TWO  NARDA  PROBES  AND  STORED  ON  DISC,  PLOT  DATA  - 

50  !  *  IN  VARIOUS  FORMS  AND  DO  DOSIMETRY  CURVE  FITTING  - 

60  !  * 

70  !  *  WRITTEN  BY  CHRISTOPHER  C  DAVIS,  DECEMBER.  1986  * 

80  !  * 

90  . . . 

100  OPTION  BASE  1 

110  MASS  STORAGE  IS  ":D500" 

120  DIM  ATEMP(  1  000  ) , ATIME< 1000) ,BTEMP( 1000  ) ,BTIME( 1000  ) , PONT  (50) ,PGFFT< 
130  DIM  ANUM(20) ,BNUM(20) , IP ON (20) ,IPOFF(20) ,  DUMMYS [ 50  3 
IaO  ♦  ***  FUNCTION  TO  GIVE  INVERSE  VIDEO  DISPLAYS 
150  DEF  FN IV IDS ( DUMMYS ) 

160  LENGTH-LEN  (DUMMYS) 

170  FOR  11=1  TO  LENGTH 

ISO  DUMMYS! II , II ] =CHRS  <NUM  ( DUMMYS [ I I . I I ]) r 1 28 ) 

190  FNIVID$=DUMMYS 
200  NEXT  I i 
2 1  0  F  N  END 
220  CLEAR 

230  DISP  FNIVIDS("DATA  DATE?" > 

240  DISP  FNl V IDS ( "DA 1 1  rORMAT") 

250  DISP  FNIVID$( "DAY  MONTH , YEAR") 

26C  IN-'L:  •  EXPDATlS 

270  DISP  FNIVIDS( "ENTER  SAMPLE  ID") 

280  INPUT  SAMPLES 
290  EXPDATAS=EXPDATES*SAMPLES 
300  ON  ERROR  GOTu  320 
8  1  0  GO  TO  3  70 

320  IF  ERRN  =67  THEN  DISP  "YOU  HAVE  ENTERED  THE  WRONG  FILE  NAME" 

330  IF  ERRN  *67  THEN  DISP  "AN  ERRjR  HAS  OCCURRlD  -  START  AGAiN" 

340  IF  ERRN  *67  THEN  GOTO  10 

350  CAT  " : D50  1  ” 

360  DISP  "These  ARE  tHE  DATA  FILES  -  CHOOSE  THE  ONE  YOU  WANT"  <  GOTO  23 
370  ASSIGN-  1  TO  EXPDATAS6" : D50 1 " 

380  NUMSMU  =  0  »  NTEST  =  0  i  SAR  =  0  »  TODS  =  "0" 

390  IhT0T  =  0  H  IBTuT  =  C  R  IONT  =  0  <s  IOFFT  =  0  £  I A 1  =  *  6  IB".  =1  C  I3N'  =  ‘  S  10” 

400  READ*  1  :  IATGT  1  READ  NO  OF  PROBE  A  READINGS 

410  FOR  I  =  I A 1  TO  IATOT  *  READ*  1  ;  ATEMP  ( I  >  ,  AT  IME  ( I  )«■  NEXT  I  !  PROBE  A 

420  READ*  1  ;  IBTOT  !  READ  NO  OF  PROBE  B  READINGS 

430  FOR  I  =  I B 1  TO  IBTOT  *  READ*  1  ;  BTEMPvI)  ,BTIME<  I)&  NEXT  I  !  PROBE  B 

440  READ*  1  ;  IONT  !  NO  OF  TIMES  MW  PWR  TURNED  ON 

450  FOR  I  *  ION  1  TO  IONT  &  READ*  1  ;  PONT(I)^  NEXT  I  !  ON  TIMES 

460  READ*  1  ;  IOFFT  !  NO  OF  TIMES  MW  PWR  TURNED  OFF 

470  FOR  I  =  IOFF 1  TO  IOFFT  9  READ*  1  ;  P0FFT(I)9  NEXT  I  !  OFF  TIMES 

480  IA1  =  IA1 +IATOT  3  IB  I  =  IB  1 +IBTOT  9  ION  1 -IQN1+IONT  &  I  OFF  1  =  IOFF 1 ♦ IOFFT 

490  ON  ERROR  GOTO  499 

491  8EAD*  1  ;  TODS 

499  Or F  ERROR  A  ASSIGN*  1  TO  * 

500  !  PURGE  EXPDATAS4" : D50 1 " 

510  DISP  "TOTAL  PROBE  A  READINGS  =  " , IATOT 

520  !  FOR  I-i  TQ  IATOT  9  DISP  ATEMP( I > , AT IME ( T ) , I  9  NEXT  I 

530  DISP  "TOTAL  PROBE  B  READINGS  -  ".IBTOT 

540  »  FOR  I-i  TO  IBTOT  9  DISP  BTEMP< I > , BTIME <  >,I  9  NEXT  I 

554  DISP  "MW  POWER  TURNED  ON  ".IONT,"  TIMES" 

S6U  3ISP  "MW  POWER  TURNED  OFF  ", IOFFT,"  TIMES" 

570  DISP  "THE  MW  PWR  WAS  TURNED  ON  AT  THE  FOLLOWING  TIMES" 

580  'OR  1=1  TO  IONT  9  DISP  PONT(I)  9  NEXT  I 
5U0  DISP  "Sample  ID  is  " ,FNIVIDS<SAMPLES > 


oOO 
601 
602 
610 
620 
630 
640 
650 
660 
670 
bSO 
£90 
700 
7  I  0 


DISP  "Date  ot'  exposure  was  " .  FNIVIDSC tXPDATES) 

IF  TOD$="0"  THEN  GOTO  610 

DISP  "Tine  of  day  was  ” . FNI V IDS< TODS > 

t  **•*«**»*•»»»***»**»»»»»**»»**»•«*»***»***»»*»**• 

!  *  SECTION  FOR  SMOOTHING 

?  «  AND  PLOTTING  OF  EXPERIMENTAL  DATA 
!  *  ■* 
!  »  WRITTEN  BY  Christopher  C.  Davis,  June  1964.  * 

!  *  MODIFIED  DECEMBER  1986 


DIM  AA ( 5 ) ,SN(21 ,21 ) , R ( 20 ) ,HEAT<5> ,CG0L(5> 


!  ARRAYS  USED  IN  CURVE  SMOCTr. 


720  NYSLS="N" 

730  IF  ATIME( I A7GT ) > BT IME ( I3TCT >  THEN  XMAX=A7  ImE ( IATQT  )  ElSE  XMAX  =  B7  IME (  1 9 T 
)  1  SET  X  AXIS  MAX  TO  MAX  OBSERVATION  TIME 
740  X M I N  =  0  1  INITIALLY  SuT  X  u  R I O I N  TO  0 

750  JM,IN=1  ti  KM  I  N  =  1 

T60  DISP  "TEMPERATURE  MEASUREMENT  INTERVAL  *AS",XMAX 

770  YMIN-10'20  *  YMAX  =  - ( 1 0 ' 20  > 

780  !  *•»  FIND  MAXIMA  AND  MINIMA  OF  DATA 
790  FOR  1=1  TO  IAT07 

300  IF  ATEMPc I )> YMAX  THEN  YMAX= ATEMP< I > 

810  IF  ATEMP<  I  X  YMIN  THEN  YMlN  =  ATEM-’<  I ) 

820  NEXT  I 

830  FOR  1=1  TO  IBTCT 

8^0  IF  STEMP<  I )  >  YMAX  Ti-iEN  YMAX  =  BTEMP  <  I ) 

850  IF  BTEMPdXYMlN  THEN  YM.lN  =  BTEMP(  I ) 

6b0  NEXT  I 

370  NUMPLT=0  4  IF  NTESTfO  THEN  CLEAR 

380  DISP  "DO  YOU  WISH  TO  PLOT  PROBE  A,  PROBE  B.  OR  BOTHA  "dFN I V I DS< "ENTER  1 
OR  3") 

390  DISP  "YOU  CAN  ONLY  CALCULATE  SAR  USING  PROBE  A" 

900  INPUT  NPROEE 

910  IF  NPRCBE= 1  THEN  GOTO  930 

920  GOTO  940 

930  FOR  1=1  TO  IATOT  4  BT IME (  I  > =AT IME<  I  >  4  BTEMP ( I )  =ATEMP< I )  4  NEXT  I 
940  IF  NTEST-0  THEN  GOTO  960 

950  NYSL3="Y" 

960  DISP  "ENTER  TIME  REGION  YOU  WISH  TO  ANALYZE" 

970  INPUT  TMIN , TMAX 

980  FOR  1=1  TO  IATOT  !  FIND  RANGE  OF  1EMP  ARRAY  TO  BE  PLOTTED 

990  IF  TMIN>ATIME( I )  THEN  JMIN=I 

1000  NEXT  I 

1010  FOR  1=1  TO  IBTQT  a  IF  TMIN>BT IME < I >  THEN  KMIN=I 
1020  NEXT  I 

1030  IF  ATIME(  JMINXBTIME(KMIN)  THEN  TMIN  =  AT  IME (  JM IN )  ELSE  TMIN  =  BTIME(KMIN) 
1040  FOR  1-1  TO  IATOT  6  IF  TMAX>AT IME ( I )  THEN  JMAX-I 
1050  NEXT  I 

1060  FOR  1-1  T-0  IBTOT  $  IF  TMAX>BTIME(  I  >  THEN  KMAX-I 
1070  NEXT  I 

1080  IF  AT IME ( JMAX  >  >BTIME(KMAX  >  THEN  TMAX  =  AT IME ( JMAX >  ELSE  TMAX  =  BT IME ( KMAX > 

1090  IF  NTEST-0  THEN  GOTO  1400 

MOO  IF  NYSL$-"N"  OR  NYSL$="NO"  THEN  GOTO  1400 

1110  !  *-*  SECTION  FOR  SELECTION  OF  CURVE  REGION  FOR  SMOOTHING 

M20  TEMAX--1  .E20  9  TEMIN-1  .E20 

1130  FOR  I  =  JM IN  TO  JMAX  &  IF  ATEMP< I XTEMIN  THEN  TEMIN-ATEMP ( I ) 

1140  iF  A  TEMP ( I ) - TEMIN  THEN  JTURN=I  ?  Find  inflection  point  on  temp  curve 
1150  NEXT  I 


I 


1160  DISP  "NUMBER  OF  POINTS  ON  COOLING  SLOPE-  " , JTURN- JMIN+ 1 

1170  DISP  "NUMBER  OF  POINTS  ON  HEATING  SLOPE-  " , JMAX- JTURN+1 

1180  DISP  "IS  THIS  SATISFACTORY?  Y(ES)  OR  N(0>"  §  INPUT  NYSATS 

1180  IF  NYSAT$-"N"  OR  NYSATS-"NO"  THEN  GOTO  960 

1200  IF  NYSLS-"N"  OR  NYSLS«"NO"  THEN  GOTO  MOO 

1210  CLEAR  «  DISP  "COOLING  REGION  POINTS" ,JMIN, JTURN 

1220  DISP  "HEATING  REGION  POINTS" , JTURN , JMAX 

1230  ON  KEY#  1, "COOLING"  GOTO  1280 

1240  ON  KEY#  2, "HEATING"  GOTO  1340 

1250  DISP  FNIVIDSC'SELECT  OPTION")  «  DISP  6  DISP 

1260  KEY  LABEL 

1270  GOTO  1270 

1280  NSM-1  £  DISP  "ENTER  POINTS  RANGE  FOR  COOLING  SLOPE" 

1290  DISP  "ENTER  0,0  FOR  HELP  IN  FINDING  TURNING  POINT  IN  HEATING  CURVE"  i  INPJT 
JMIN, JTURN 

1300  IF  JMIN=0  AND  JTURN-0  THEN  GOSUB  Find 
1310  DISP  "JTURN-", JTURN 

1320  DISP  "ENTER  JMIN. JTURN"  £  INPUT  JhIN, JTURN 
1330  GOTO  1390 

1340  NSM  =  2  £  DISP  "ENTER  POINTS  RANGE  FOR  HEATING  SLOPE" 

1350  DISP  "ENTER  0,0  FOR  HELP  IN  FINDING  TURNING  POINT  IN  HEATING  CURVE"  5  IN=»J' 
JTURN,  JMAX 

1360  IF  JTURN-0  AND  JMAX=0  THEN  GOSUB  Find 
1370  DISP  "JTURN-", JTURN 

1380  DISP  "ENTER  JTURN. JMAX"  £  INPUT  JTURN, JMAX 

1390  NUMPlT=1  6  GOSUB  Smooth 

1400  GOSUB  Aplot 

1410  ON  KEY#  1, "AGAIN"  GOTO  1530 

1420  ON  KEY#  2, "DOSE"  GOTO  1580 

1430  ON  KEY#  3, "PLOT"  GOTO  1590 

1440  ON  KEY#  4 , "SCRDMP"  GOTO  3910 

1450  ON  KEY#  5. "EXIT"  GOTO  3920 

1460  CLEAR  *  KEY  LABEL 

1470  DISP  "SELECT  OPTION  -  "aFNlVlD$("DOSE" >&"  ALLOWS  CALCULATION  OF  SAR" 

1480  DISP  "  -  "&FN IV IDS ("AGAIN" )&"  RE-READS  DATA  FRQM  DISC" 

1490  DISP  "  -  "&FN IVIDSC "PLOT” )  &  "  RE-PLOTS  LAST  PLOT  -  USE  WITH  C 

ARE" 

1500  DISP  "  -  "&rN IV IDS ("SCRDMP") a"  DUMPS  GRAPHICS  TO  PRINTER" 

1510  DISP  "  -  "&FNIVIDS("EXIT")&"  TERMINATES  EXECUTION- 

1520  GOTO  1520 

1530  CLEAR 

1540  DISP  "IF  YOU  CONTINUE  DATA  WILL  BE  RE-READ  FROM  DISC" 

1550  DISP  "PRESS  "4FNIVIDS( "AGAIN" )a"KEY  IF  YOU  WISH  TO  DO  THIS" 

1560  ON  KEY#  1, "AGAIN"  GOTO  370  £  KEY  LABEL 
1570  GOTO  1570 
1580  GOTO  930 

1590  IF  NUMSMO-O  THEN  GOTO  1620 

1600  GOSUB  Aplot  £  GOTO  1520 

1610  GOTO  1410 

1620  DISP  "OPTION  NOT  AVAILABLE"  £  GOTO  370 
1630  Smooth: 

1640  f  . . ****** 

1650  !  •  • 

1660  !  •  CURVE  SMOOTHING  SECTION 

1670  f  *  * 

1630  ! 

1690  DISP  "CURVE  FITTING  IN  PROGRESS  "6FNIVIDS("PLEASE  WAIT") 

1700  NUMSMO-1 

1710  IF  NSM-2  THEN  GOTO  1750 

1720  N* JTURN- JMIN+ 1 

1  730  FOR  I  -  !  TO  N  £  BTIME(  I ) -BT  IMF  ( JMTN  +  I  -  1  )  £  BTEMPf  I )  ••BTEMP<  JMIN  +  I  -  1  >  £  NE  * 


T  I 
I  7m0 
•  75u 
1  760 
RN  + 1 -  1  ) 
1770 
1  760 
!  790 
1300 
1610 
1620 
1830 
1890 
1850 
1  dbO 
1670 
:  660 
;  DbO 
’  9  1 0 
'  9 ' 0 


.i  u 

‘.-50 


N  , 


GOTO  1780 
N=JMAX~ JTURN+1 

F3R  1  =  1  TO  UMAX- JTURN+  1  C*  B  T I  ME  <  I >  =  BT IME ( JTURN  + 1 -  1  )  #  BTEMP<  I )  =BTEMP( JTL 


NEXT  I 

FOR  1=1  TO  5 

R(I)-0  *  AA<I>=0  »  FOR  J=1  TO  21  j»  SN(I,J>=0  @  NEXT  J  ■*  NEXT  I 

DISP  "WHAT  ORDER  OF  POLYNOMIAL  DO  YOU  WISH  TO  FIT  -  1-LINEAR.. 

INPUT  NS 

FOR  J-1  TO  NS  +  1 

FOR  1=1  TO  N 

IF  J  =  NS  + 1  AND  B T I ME < I »  =  0  THEN  GOTO  1860 
R(J)=BTEM^< I)»BTIME(I) ' <NS+1-J)+R< J)  6  GOTO  1870 
R< J)=BTEMP( I ) +R ( J  ) 


NEXT  I 
NEXT  J 

FOR  K  = 1  TO  NS  + 1  »  FOR  J=1  TO  NS+1  ?  FOR  1=1  TO  N 
n  2«NS-J-N+2=0  AND  BT IME< I ) =0  THEN  GOTO  1920 
SN(o ,n)  =  B  T  I  ME  <  I  )  ‘  <2*NS-J-<+2>+SN<J.K>  -  GOTO  1920 

sn<w  .»:>  =  ’  +sn«.  j .  k  > 

NEXT  I 
NE>.  T  j 
NEXT  < 


«« 


rwK  1=1  TO  NS+1  3»  FOR  J=l  TO  NS+1 

AH  (  N;  ♦  _  -  I  )  fc$N  <  NS  +  2-  I  ,  J  )  +  R  (  J  >  +  AH  (  NS  +  2  -  I  >  A  NEXT  J  *  f,E>”  I 

IF  NSM=1  ThEN  DX=ATIME<  jTuRN  )  -  AT  IME  «.  JMIN)  ELSE  DX  =  A  T  IME  (  UMAX  ) -A  T I 0 


.A  „  =  1 


J  0  5  0 
6  6  6  0 
6  T  7  j 
-  -J 


- :  oa 
2:10 
2120 
2130 
:  i^o 
2150 
2160 
2170 
2180 
2190 
2200 
2210 
2220 
2230 
2290 
2250 
2260 
2  27  0 
22b0 
2290 

2  ;f  "j 


IF  NGM-1  THEN  X$  T  ART  =  AT I  ME  <  JMI N )  ELSE  XST  RT-ni  iflE(JTuRN) 

+  06  N=i  TO  2GC  I  =  K  +  800  t-  BTlMidl  w.-l  )*DX/20T)  +  xSTART  «  STEM-"  ( 

TO  NS+1 

IF  NS+1-J+0  AND  5  T I  ME  <  I )  =  0  TnE‘i  GOTO  2050 
BTE*P<  I )  =  A  A  (  J  )  *EI  IME(  I  >  (NS+  t  -  j  )  +B  :lN-  ,  *  '■ 
o  u  1  0  2  0 1.  J 

BTEM^( I )  =  AA( J ) +B'£HP(  I  > 

NEXT  J 
NlX*  + 

RE  T  ^,n 

I  ri  ■/ 1  : 


MATRIX  .  N  V  l  A  3 1  u  N  BY  GAUSS  -  JCRD  AN  EL  J.MINAT  ION  WITH  FULL 
SQUARE  MATRIX  Ur  ORDER  M~ 


PILOTING 


D .  :1  rP(5),oP<5) 

NP-NS+ 1 
NCQL=NP 

INITIALIZE  PERMUTATION  LECTORS 

FOR  1=1  TO  MP  #  KP ( I )  =  I  9  JPU)  =  I  *  NEXT  I 

INVERT  MATRIX 

FOR  IR-!  TO  NP 

FIND  THE  PIVOT  ELEMENT 

PlyoT-o 

FOR  I  =  I R  TO  NP  «  FOR  J  =  I R  TO  NP  >  K  =  KP<I)  •*  L  =  JP(J) 
TEST  =  ABS  (SNO  ,L)  )-AB5  (PIVOT) 
i  T  Tr  5  T  <  -  ij  tmEN  GO  1  0  23  7  0 


)  =  'j  6  : 


2320  IF  TEST >0  THEN  GOTO  2340 

2330  GOTO  2370 

2340  IK-I 

2350  JK  =  J 

23S0  PIVOT-SMK.L) 

2370  NEXT  J  6  NEXT  I 

2380  ! 

2380  !  THE  PIVOT  ELEMENT  IS  SN(K,L) 

2400  ! 

2410  !  UPDATE  PIVOT  VECTORS 

2420  ! 

2430 
L 

2^40  ! 

2450  ! 

2480  ! 

24  70 

2  4  S  0  ! 

2480  ! 

2500  * 

25  >0 
2520 
253U 
2540 
2550 
j5o  0 
25  70 
2530 
2530 
2800  ! 

2610  ! 

2820  ! 

“t  r  v .  • 

u 

2'  b  5 

.vjOU  b. 

-570 

4  d  u  :  *  » 

-c-30  !  *  PLOTTING  SUBROUTINE  * 

2 /  —  *  —  —  —  *  —  —  —  —  —  * 

2  7 1  0  ! 

2720  NlT=1 

2730  CLEAR  i<  IF  NUMPLT=0  THEN  GOTO  2930 
2740  ALPHA  *  DISP  "PARAMETERS  OF  FIT"  £  DISP 
2750  t-GR  1=1  TO  NP  «t  DISP  AA(I)  «  NEXT  I 
2780  IF  NSM  =  2  THEN  GOTO  2790 
2770  FOR  1  =  1  TO  NP  «  COOL  ( I )  =AA(  I )  <«<  NEXT  I 
2780  GOTO  2800 

2790  FOR  1  =  1  TO  NP  $  HEAT < I ) =  AA( I )  $  NEXT  I 
2800  CHI2-0 

2810  IF  NTEST<2  THEN  GOTO  2850 
2820  SAR-4 .2*(HEAT < 1 )-COOL(1 >>*1000 

2830  IF  NS=2  THEN  SAR=2*SAR«AT IME < JTURN) +4 . 2* ( HEAT (2 ) -COOL < 2 > ) 
2840  DISP  FNIVID$("SAR  IS" ) , SAR , "W/k g" 

4850  FOR  1-1  TO  N  e  Y2=0  *  FOR  J-1  TO  NS+J 
2860  IF  NS  + 1  - J-0  AND  BTIME<I)-C  THEN  GOTO  2880 
2o/0  Y2  =  AA<  J;«BTII1E<  I )  <NS+1-J)*Y2  »  GOTO  '’390 
2980  Y2-AA<  J)  +  Y2 
2890  NEXT  J 

4S-U0  CHI2=(BTFMP(  I )  -  Y  2 )  '2/Y2  +  CHI2 
29-  0  N-xT  I 


K*KP ( IK  >  ■»  <P(  IK )  =KP(  IR)  KP(IR)=K  A  L  =  JP(jK)  »  JP  ( JK  )  =  JP  ( IR )  A  JP< 


COMPUTE  NEW  ELEMENTS  OF  PIVOTAL  ROU 

FOR  J=1  TO  NCOL  A  SN( K , J ) =SN( K , J > /P I VOT  A  NEXT  J 

COMPUTE  REMAINING  ELEMENTS  OF  IR'TH  STEP 

FOR  1=1  Tu  NP 

IP  I  -  K  =  0  THEN  GOTO:  2590 

A i!_  =  5N (  I  ,i_) 

SN( I  ,L  )  =  - caIl/p  IVCT  ) 

FOR  J=1  TO  NCOL 

IF  j  * l  =  0  THEN  G-TC  2580 

SN(  I ,  j  )  =  SN(  I  ,  J  )  -  AIL* Si.  (.<  ,J) 

NEXT  j 
NEXT  I 


INVER7  PIVOT  E-EMENT 

Sit  <  K  .  l  )  =  !  /  P  I  VC ”  «  NEXT  IR 
RE  7  URN 


2920  DISP  “OBSERVED  CHI2=  ’*  .CHI2 
2930  DISP  "DO  YOU  WANT  A  HARD  COPY 


WITH  "  ,  M  +  NS+  1  D.F." 

PLOT'?  Y  < ES )  OR  N ( 0 ) " 


2540  INPUT  NYSi*  IF  NUMSMO=1  THEN  GOTO  3040 

2950  DISP  "TIME  KILL  BE  PLOTTED  FROM  " , XMIN , "TO" , XMAX 

2960  DISP  "TEMPERATURE  WILL  BE  PLOTTED  FROM"  ,  YMIN  , "TO** ,  YMAX 

2970  DISP  "IF  YOU  ARE  SATISFIED  PRESS  RETURN" 

2380  INPUT  NYLIMS 

2990  IF  NYLIMS*""  THEN  GOTO  3040 

3000  DISP  "ENTER  TIME  RANGE  TO  BE  PLOTTED,  TMIN.TMAX" 

3010  INPUT  XMIN, XMAX 

3020  DISP  "ENTER  TEMPERATURE  RANGE  TO  BE  PLOTTED" 

3030  INPUT  YMIN, YMAX 

3040  IF  NTEST>=  1  THEN  GOTO  3340 

3050  CLEAR 

3060  DISP  "CALCULATING  PLOT  SCALING  PARAMETERS  "AFN IV  IDS ( "PLEASE  WAIT") 
3U70  AXMAX= INT  ( LGT  (XMAX)) 


3080  AYMAX=INT  (LGT  (YMAX)) 

3090  XF=  1  0  ’  AXMAX  «,  YF=  1  0  'AYMAX  £  DELX-. 01  «  DELY-.G1 
3100  XPWRS-VALS  (AXMAX)  *  YPWRS-VAL3  (A YMAX ) 


3110  !  XMIN-0  s  Y M I N  =  0 


3120 

IF 

( XMAX-XMlN  > /XF>  = 

.  1 

THEN 

QELX=.G2 

3130 

IF 

( XMAX-XMIN ) /XF>  = 

.3 

THEN 

DELX*  .05 

3140 

IF 

( XMAX-XMIN ) /XF>  = 

.5 

THEN 

DELX= . 1 

3150 

IF 

( XMAX-XMIN) /XF>= 

1 

T  HEN 

DlLK=  .2 

3180 

IF 

(XMAX-XMIN)/XF>= 

3 

THEN 

DELX=  .  5 

31  70 

IF 

(XMAX-XMIN )/XF>= 

5 

i  HEN 

D  E  L  X  =  1 

3180 

IF 

( YMAX -YMlN)/YF>= 

.  1 

THEN 

QELY- .02 

3190 

IF 

( YMAX- YM IN ) / YF>= 

.3 

THEN 

DElY=  .05 

3200 

IF 

( YMAX- YM I N ) / YF  >  = 

.5 

THEN 

DEl  Y  =  .  1 

3210 

IF 

( YMAX-YMIN)/YF>- 

1 

T  HEN 

DELY=  .2 

3220 

IF 

( YMAX- YMIN) /YF  >  * 

3 

THEN 

DEL  Y  =  .  5 

3230 

IF 

<  YMAX-YMIN) /YF>  = 

5 

THEN 

DEL  Y=  1 

32^0  DISP  "THE  X-AXIS  TICK  SPACING  WILL  BE” , DELX 

3250  DISP  "IF  YOU  ARE  SATISFIED  PRESS"&FN I VIDS( "ENDlINE"  ) 

3260  INPUT  NIKS*  IF  NIKS-""  THEN  GOTO  3280 

3270  DISP  "ENTER  X-AXIS  TICK  SPACING"  *  INPUT  DElX 

3280  DISP  "THE  Y-AXIS  TICK  SPACING  WILL  BE" ,  DEL Y 

3290  DISP  "IF  YOU  ARE  SATISFIED  PRESS"£FNI VIDS ( “ENDLINE" > 

3300  INPUT  NTKS£  IF  NTKS-""  THEN  GOTO  3320 

3310  DISP  "ENTER  Y-AXIS  TICK  SPACING"  *  INPUT  DElY 

3320  XNUM  =  CEIL  ( (XMAX-XMIN)/(XF*DELX> ) 

3330  YNUM  =  CEIL  <( YMAX- YMIN )/( YF-DELY > ) 

3340  IF  NY$="N"  THEN  GOTO  3330 
3350  PLOTTER  IS  505 

3360  GRAPH  $  IF  NUMPLT-0  THEN  GCLEAR  *  DEG 
3370  GOTO  3420 
3380  PLOTTER  IS  1 
3390  GRAPH 

3400  !  IF  NUMPLT*0  AND  NTEST>0  THEN  GOTO  3110 

3410  IF  NUMPLT-0  THEN  GCLEAR 

3420  LOCATE  30, RATIO  *100-10,20,95 

3430  !  IF  NUMPLT#D  AND  NTEST>0  THEN  GOTO  3110 

3440  SCALE  XMIN/XF , XMIN/XF* XNUM*DELX , YMIN/YF , YMIN/YF+ YNUM*DELY 

3450  IF  NTEST>0  AND  NUMPLT*0  THEN  GOTO  3870 

3460  FXD  2 

3470  LGRID  DELX , DEL Y , XMIN/XF , YMIN/ YF 
3480  ftXS=(XMAX-XMIN)/(XF*1 0) 

3490  AYS* ( YMAX- YMIN) / ( YF* 1 0 ) 

3500  MOVE  (XMAX*XMIN)/<2*XF)  ,  YMIN/ YF -2-AYS  *  LORE.  6 
3510  LABEL  USING  "K"  ;  "TIME  (SECS)  x  lE"AXPWf<5 
3j/0  MOVE  XMIN/XF-2-AXS, ( YMAX+ YM IN ) / ( 2- YF ) 


3530  LQRG  4  $  DEG  3  LDIR  90 

3540  LABEL  USING  ,,K"  ;  "TEMP .  (DEGS  C)  x  1E"&YPWR$ 

3550  LDIR  0 

3560  IF  NPR0BE=2  THEN  GOTO  3630 
3570  LINE  TYPE  1 

3580  MOVE  ATIME( JMIN>/XF , ATEMP< JMIN) /YF 
3590  FOR  I-JMIN  TO  JHAX 
3600  DRAW  ATIME<I)/XF,ATEMP(I)/YF 
3610  NEXT  I 

3620  IF  NPROBE* 1  THEN  GOTO  3660 
3630  LINE  TYPE  5 

3640  MOVE  BTIME(KMIN) /XF . BTEMP (KM IN)/ YF 

3650  FOR  I-KMIN  TO  KMAX  $  DRAW  BT IME < I > /XF , BTEMP < I ) / YF  «  NEXT  I 
3660  IF  NYSL$="N"  OR  NYSLS="NO"  THEN  GOTO  3720 
3670  LINE  TYPE  4 

3680  MOVE  BT IME  <  80 1 ) /XF , BTEMP ( 80 1 ) /YF 
3690  FOR  1=801  TO  1000 
3700  DRAW  BTIME<I)/XF.BTEMP(I)/YF 
3710  NEXT  I 


3720  PEN  UP 

3730  MOVE  XMIN/XF , YMIN/ YF 

3740  IF  SAR  =  0  THEN  GOTO  3780 

3750  MOVE  XMIN/XF  +  4 . 5*AXS  ,  YMAX/YF+.NAYS 

2760  LORG  4 

3770  LABEL  "SAR", SAR , "W/k g" 

3760  NTEST=NTEST +  1 

3735  DISP  "PRESS  “iFNIVIDS("CCNTINUE">£"TQ  CALCULATE 
3790  PAUSE 
3800  KEY  LABEL 
3810  RETURN 


38c 0  Find:  *  ***  SUBROUTINE  TO  FIND  POINT  WHERE  HEATING  STARTS 
3330  D.SP  "ENTER  A  T E S T  RANGE  J1,J2"  4  INPUT  J1.J2 
3840  DISP  "TIME"  TEMP"  ,  "DT/'Dt " 

3850  FOR  I=J1  TO  J2-‘ 

3860  DISP  ATI  ME (I) , ATEMP< I > , ( ATEMP ( I  +  1 ) - ATEMP ( I ) ) / ( AT IME ( I  +  1 ) -AT IME < I ) ) 
ob70  NlXt  I 


disp  "try  another  range  -  r(ES>  or  n<o>"  «  input  nyrs 

3680  IF  NYPS="Y"  OR  NYRS=" YES"  THEN  GOTO  3830 
3900  DISP  "ENTER  JTURN"  i  INPUT  jTuRN'**  RETURN 

3910  Chain  "dump" 

2920  END 


PROGRAM  FOR  DATA  ANALYSIS,  CURVE  SMOOTHING  AND  GRAPHICS  -  version  B 


U 

0 

0 

0 

0 

0 


•  *■ 

-  PROGRAM  TG  READ  TEMPERATURE  READINGS  TAKEN  FROM  * 

*  TWO  NARDA  PROBES  AND  STORED  Ol.  DISC.  PLOT  DATA  * 

*  IN  VARIOUS  FORMS  AND  DO  DOSIMETRY  CURVE  FITTING  * 

*  f* 

*  WRITTEN  BY  CHRISTOPHER  C  DAVIS.  DECEMBER.  1936  * 


0  0 
ID 
CD 
30 
-JO 
50 
id 
70 
SO 
50 
00 
1  0 
:o 

30 

AO 

50 

6  o 

70 

80 

Qj 

uO 
’0 
2  0 
3'J 
"|J 
50 

C’U 

70 

BO 

90 

GO 

*  0 

20 
30 
40 
TO 
,0 
70 
TO 
40 
JO 
10 
10 
i  0 
40 

j') 
0 
'  u* 

'/J 


" : D500" 

_  .  A  TIME ( 1000  ) .9 TEMP ( 1000  ) .BTIME< 1 000  )  ,PGNT( 50) 

ANUM  ( 20 )  .  BNUM  <  20 )  ,IPON(20>  ,IPGFF(20)  .BUM*YS[50I 

INVERSE  VIDEO  DISPLAYS 


OPTION  BASE  1 
MASS  STORAGE  IS 
DIM  ATEMP( 1 000) 

DIM 

?  ***  function  to  give 

DEE  FNIVIDS ( DUMMYS  > 
LENGTH=lEN  (DUMMYS) 

FOR  11=1  TO  LENGTH 
DUMMYS[ 1 1 . 1 : ] =  CHKS 
FN I VIDS “DUMMYS 
NEXT  II 
FN  END 
CLEAR 

DISP  FNIVIDSC "DATA 
DIS?  FNI V IDS ( ”DA~E 
DISP  FNIVIDS( "DAY 
INPUT  EXPDATES 
DISP  FNI V I DS  (  " ENTER  SAMPLE 
INPUT  $AM “’LL S 
EXPDA T AS “EXPDATESdSAMPLES 


'LEFT  (  5 J > 


(NUM  ( DUMMY'S!  I 


1 3 )  1  2  S ) 


Da  ' 


) 


c.  ' 

FORMAT”  ) 
MGNT ri ,  YEAR”) 


ID") 


ASSIGN*-  1 
NUMSMQ=0 
I A  T  G  T  =  0  s» 
RuAD*-  1  ; 

F  jR  I  =  i  A.  1 
READ*-  1  ; 
FOR  1=131 
READ*  1  ; 
FuR  I  =  1 0 N I 
READ*-  1  : 


TO  EXPDATASa” : DSC i 
v  NTEST  =  0  i  SAR  =  0 
I3TGT  =  0  *  I G N T  =  0  3 
IATuT  !  READ  Nu  OF 
TG  IATOT  i  READ*  i 
IBTuT  !  READ  NO  Or 
TG  I3TGT  i>  READ*  1 
IGN7  !  NG  OF  T IMES 
TO  IGNT  i  READ*  1 
IOFFT  *  NO  OF  TIMES 


:  b  i  =  t 


NEX' 


IOFFT  «  READ* 

4  I B 1 = I B 1 +IBTOT 


IGFF  T  =  U  j»  IA1  =  1 
PROBE  A  READINGS 
:  A  TEMP ( I ) , A  T  I  ME ( I )  ? 

PRUBt.  a  READINGS 
;  BTEMP(I) , B T IME ( I ) 3  NEXT  I  !  P 
Mh  PWn  TURNED  ON 
;  PONTOi*  NEXT  I  !  ON  TIMES 
MW  PWP  TURNED  OFF 
1  ;  PO)FT<I>»>  NEXT  I  !  OFF  TIME 
(j  ION  1  =  ION'  +  IONT  G  IOFFT  =  IDFF  1 


IATOT 


DISP  ATEMP(I) ,ATIME(I).I  «  NEXT 


FOR  I-IGFF1  TO 
IA 1 =  I A 1 + 1 ATOT 
ASSIGN*  1  TO 
!  PURGE  EXPDATAS6" :D501 " 

DISP  "TOTAL  PROBE  A  READINGS 
!  FDR  1=1  TO  IATOT  G 
DISP  "TOTAL  PROBE  B  READINGS  -  " , IBTOT 
!  FOR  1-1  TO  IBTOT  G  DISP  B TEMP ( I )  .  BT IME ( I ) 

DISP  "MW  POWER  TURNED  UN  ",IONT,"  TIMES" 

DISP  "MW  POWER  TURNED  OFF  ".IOFFT,"  TIMES" 

DISP  "THE  MW  PWR  WAS  TURNED  ON  AT  THE  FOLLOWING  TIMES 
FOR  1-1  TO  IONT  G  DISP  PONT(I)  G  NEXT  I 
DISP  “Sample  ID  is  " ,FNIVID$(SAMPLES) 

DISP  "Date  of  exposure  uas  ", FNIVID$< EXPDATES  ) 


ION  1  =  1  i  IGFF  *  = ! 


h ' j  £  c_  A  DA  -  A 
Rur E  3  DA  A 


.  ur  r  ! 


I  G  NEXT  I 


SECTION  FOR  SMOOTHING 
AND  PLOTTING  OF  EXPERIMENTAL  DATA 


iRITTEN  BY  Christopher  C.  Dqvis 
MODIFIED  DECEMBER  I'bnF. 


June  1984 


6io  ! . * . * . -*  -**•*.* . . — 

620  ! 

630  ! 

6^*0  DIM  AA<5> ,SN(21 ,21 > ,R(20>  ,  HEAT <5>  ,C00L(5>  !  ARRAYS  USED  IN  CURVE  SMOGTH! 

650  NYSLS-"N" 

660  IF  ATIME(IAT0T>>BTIME(IBT0T)  THEN  XMAX-ATIME< IATQT)  ELSE  XMAX“BTIMM.C <  IBTC 
)  !  SET  X  AXIS  MAX  TO  MAX  OBSERVATION  TIME 
670  XMIN-0  !  INITIALLY  SET  X  ORIGIN  TO  0 
680  JMIN-1  6  KMIN-1 

690  DISP  "TEMPERATURE  MEASUREMENT  INTERVAL  WAS",XMAX 
700  YMlN-10‘20  G  YMAX-- < 1 0 '20 > 

710  !  •**  FIND  MAXIMA  AND  MINIMA  OF  DATA 
720  FOR  1-1  TO  IATOT 

730  IF  ATEMP( I)>YMAX  THEN  YMAX-ATEMPC I > 

740  IF  ATEMP< I ) <YMIN  THEN  YMIN*ATEMP( I ) 

750  NEXT  I 

760  FOR  1-1  TO  IBTOT 

770  IF  BTEMP(I)>YMAX  THEN  YhAX-BTEMP< I ) 

780  IF  BTEMP< I  X YMIN  THEN  YMIN-BTEMP< I ) 

780  NEXT  I 

600  NUMPLT-0  3  IF  NTEST#0  THEN  CLEAR 

810  DISP  "DO  YOU  WISH  TO  PLOT  PROBE  A,  PROBE  B,  OR  BOTH?  "&FNIVIDS("EN  ~EP  1, 
OR  3") 

820  DISP  "YOU  CAN  ONLY  CALCULATE  SAR  USING  PROBE  A" 

830  INPUT  NPROBE 

840  IF  NPROBE* 1  THEN  GOTO  860 

850  GOTO  870 

860  FOR  1-1  TO  IATOT  *  BTIME<  I ) -HTIMEU )  *  BTEMP(I>-ATEMP(I>  «  NEXT  I 
870  IF  NTEST-0  THEN  GOTO  890 

880  NYSL$-"Y" 

890  DISP  "ENTER  TIME  REGION  YOU  WISH  TO  ANALYZE" 

S00  INPUT  TMIN.TMAX 

910  FOR  1-1  TO  IATOT  !  FIND  RANGE  OF  TEMP  ARRAY  TO  BE  PLOTTED 
920  IF  TMIN>ATIME(I)  THEN  JMIN-I 

930  NEXT  I 

940  FOR  1-1  TO  IBTOT  d  IF  TMIN>BTIM£(I>  THEN  KMIN-3 

950  NEXT  I 


960 

970 

980 


IF  ATIME ( JMINXBTIME(KMIN)  THEN  TM]N-ATIME(  JMIN)  ELSE  TMIN=BTIME<,  <M 
FOR  1-1  TO  IATOT  e  IF  TMAX>ATIME(I>  THEN  JMAX-I 
NEXT  I 


990 
1000 
1010 
1020 
1030 
1040 
1050 
1060 
1070 
1080 
1090 
1100 
1110 
1120 
1130 
1140 
1  !  50 
1160 
1170 
1180 


FOR  1-1  TO  IBTOT  e  IF  TMAX>BTIME< I)  THEN  KMAX-I 

IF  ATIMEC JMAX > >BTIME<KMAX >  THEN  TMAX-ATIMEl JMAX)  ELSE  TMAX-3TIME  ;KMAX> 

IF  NTEST-0  THEN  GOTO  1250 

IF  NYSL$-"N"  OR  NYSLS-MNO"  THEN  GOTO  1250 

SECTION  FOR  SELECTION  OF  CURVE  REGION  FOR  SMOOTHING 
TEMAX--1 .E20  »  TEMIN-1 .E20 

ESR*fcJMlN  T0  J^AX  G  IF  ATEMP< I XTEMIN  THEN  TEMIN-ATEMP( I ) 
NEXTTIMP(I>‘TEMIN  JTURN-I  !  Find  inflection  point  on  temp  carve 

DISP  "NUMBER  OF  POINTS  ON  COOLING  SLOPE-  " , JTURN-JMIN-M 

DISP  "NUMBER  OF  POINTS  ON  HEATING  SLOPE-  ", Jt AX-JTURN+1 

DISP  "IS  THIS  SATISFACTORY?  Y<ES>  OR  N(0>"  G  NPUT  NYSATS 

IF  NYSAT$-"N"  OR  NYSAT$-"NO"  THEN  GOTO  B90 

IF  NYSL$-"N"  OR  NYSL$-"NO"  THEN  GOTO  1250 

CLEAR  &  DISP  "COOLING  REGION  POINTS" , JMIN, JTURN 

DISP  "HEATING  REGION  POINTS" , JTURN, UMAX 

ON  KEY#  1, "COOLING"  GOTO  1210 

ON  KEY#  2, "HEATING"  GOTO  1230 

DISP  FNJVIDS<"SELECT  OPTION")  6  DISP  e  DISP 


1200  GDTu  1200 

1210  NSM-I  #  DISP  "ENTER  POINTS  RANGE  FOR  COOLING  SLOPE"  «  DISP  £  DISP  »  INPU 
T  JMIN.JTURN 
1220  GOTO  1240 

1230  NSM-2  I*  DISP  "tNTER  POINTS  RANGE  FOR  HEATING  SLOPE"  y  DISP  9  DISP  I<  INPU 

T  JTURN ,  JMAX 

1 240  NUMPLT* 1  $  GOSUB  Smootn 
1250  GOSUB  Ap lot 
1260  ON  KEY#  1, "AGAIN"  GOTO  1380 
1270  ON  KEY#  2. "DOSE"  GOTO  1430 

1280  ON  KEY#  3, "PLOT"  GOTO  1440 

1290  ON  KEY#  4,"SCRDMP"  GOTO  3530 
1300  ON  KEY#  5. "EXIT"  GOTO  3540 
1310  CLEAR  «  KEY  LABEL 

1320  DISP  "SELECT  OPTION  -  "&FNIVIDS< "DOSE" M"  ALLOWS  CALCULATION  OF  SAR" 

1330  DISP  "  -  "&FNI VIDS ( "AGAIN" )6"  RE-READS  DATA  FROM  DISC" 

1340  DISP  "  -  "aFNlUIDS<  "PLOT"  )<i"  RE-PLOTS  LAST  PLOT  -  USE  Wl^H  C 

ARE“ 


’350  DISP  "  -  " AFNI V I DS( "SCRDMP" ) A"  DUMPS  GRAPHICS 

I3o0  DISP  "  -  "dFNl V  IDS ( "EXI T" )&"  TERMINATES  EXECl 

1  370  Gu  i  0  1370 
1  380  CLEAR 

1390  DISP  "IF  YOU  CONTINUE  DATA  HILL  BE  RE-READ  FROM  DISC" 

1400  DISP  "PRESS  "6r Nl VI D£< "AGAIN" )&"KEY  IF  YOU  NISn  TO  DO  THIS" 
1410  ON  KEY#  1, "AGAIN"  GOTO  300  i>  KEY  LABEl 
1420  GOTO  1420 
1  430  Gu  1 0  86  0 


0  P  R I N 


’  R" 


14^0  IF  NUMS,MO  =  0  TnEN  GOTO  1470 
1450  GOSUB  Apiot  GOTO  13  7  0 
1460  GOTO  1260 

I47j  DISP  "OPTION  Nur  AVAILABLE"  *  GOTO  300 
14o7  Snoot n : 

1  4b0  !  . . . . 

1  S  G  b  1  • 


lull. 

1520 
1530 
1540 
i  550 

1 5o0 

1570 
1  560 
T  I 
1  590 
1 60  C 
1610 


CURVE  SMOOTHING  SECTION 


DISP  "CURVE  FITTING  IN  PROGRESS  "&FN I V I DS ( "PlEASE  hAIt"> 

NUMSMO* 1 

IF  NSM-2  THEN  GuTu  1600 
N= JTURN- JMIN  + 1 

FOR  1*1  TO  N  »  BTI  ME  <1)^811  ME  (JMIN  +  I-1  )  i>  BTEMP  ( I )  =  BT  EMP<  JllIN*-]  -  1  )  e*  HEX 

GOTO  1 6  J 0 
N* JMAX- JTURN+ 1 

FOR  1*1  TO  JMAX- JTURN+  1  BT IME  <  I  >  *BT  IME  ( J  TURN- I  -  1  )  4  B  TEMP  ( I  )■  BTEMP  (  J  TU 


1620 
1630 
1640 
1650 
1660 
1670 
1680 
1690 
1700 
1  71Q 
1  720 
1  730 
1  74(J 
1  /Si 
1  767; 


C 

-J 


NEXT  I 

FOR  1-1  TO  _ 

EtU'P.  -  AA<I>"°  *  3=1  TO  21  3  SNC I ,  J)  =  0  «  NEXT  J  §.  NEXT  I 

DISP  WHAT  ORDER  OF  POLYNOMIAL  DO  YOU  WISH  TO  FIT  -  1-lINEak.. 
INPUT  NS  . 

FOR  J-1  TO  NS+ 1 
FOR  1*1  TO  n 

IF  J“NS+ 1  AND  BTIME<I)=0  THEN  GOTO  1710 
R<J>-8TEMP(1>*BTIME<I>' <NS+1 -J)+R( J)  «.  GOTO  1720 
P<J)»BTEMP< I > ♦«< J) 

NtXT  I 
NEXT  J 


FOR  r<  *  1 
IF  2*NS 


N  i 


TO  NS>1  a  FOR  J-1  TO  N5-*  !  •>  FOR  1=1 
J-K+2-0  AND  BTIME(I>-U  IntN  GLmw  i  / 


i  '<>  T  ’  '•(  <  j  >  (  '->«  v,  j 


i  U 
'u 


N 


1  > 


N<  1 


»?7C 
I  780 
1790 
1800 
1810 
1820 
1830 
1840 
N) 

1850 
1860 
DR  J« 
I  370 
i  b«U 
1390 
1900 
1910 
1920 
I  920 
1  9nC 
1950 
1360 
1970 
i960 

1  990 
oGGG 
2010 
2020 
2030 

2  U  4  C 
2050 

OOtw 

20  70 
2030 
2  090 

21  00 
2  110 
2  1  2  0 
2  1  30 
2140 

^  i  DU 
2  1  60 
2170 
2130 
2190 
2200 
2210 
2220 
2230 
2240 
2250 
2260 
2270 
2280 

2290 
2  OoO 

23  1  U 
2220 
2330 
2  d  H  0 


SNiJ.K)-  WSmU.K) 

NEXT  I 
NEXT  J 
NEXT  K 
GOSUB  In v t 

FOR  I '  I  TO  NS  + 1  «  FOR  J-1  TO  NSH 

AA<NS+2-I>“SN(NS+2-I,J)*R<J>+AA<NS+2-I>  &  NEXT  J  £  NEXT  I 

IF  NSH-1  THEN  DX*ATIME(JTURN>-ATItiE<  JttlN)  ELSE  DX=AT  IMEC  JMAX ) -AT  IME  t  JU 

IF  NSM  =  1  THEN  XSTART*ATIME<JMIN>  ELSE  XS7ART=ATIME( JTURN) 

FOR  K  =  1  TO  200  6  I-K  +  60G  $  BT  IME  <  I )  -  <  K- 1  )  *DX/20G+XS7  ART  6T£MPtI>«U  t 
I  TO  NS+ 1 

IF  NS  + 1  - J  =  0  AND  BT  IME ( I ) *0  THEN  GOTO  1900 
BTEMP(  I)*AA(J)*bTIME<l)'  ( NS  + 1  - j > +b i  tflPC  i  ) 

GOTO  1910 

BTEMP(I)=AA( J)+BTEMP( I) 

NEXT  J 
NEXT  K 
RETuRN 
Invt : 

•  MATRIX  INVERSION  3Y  GAUSS- JORDAN  ELIMINATION  HlTh  FoLl  PIvlTINu 
!  SQjhRc.  Mh  i  r\iX  Or  ORDER  NP 

t 

DIM  K  P<  5 1 , JP(5) 

NP -  N$  + 1 
NCuL  =  NP 

f 

!  INITIALIZE  PERMUTATION  sectors 

* 

FOR  I-l  Tu  NP  j  \PO*I  j-  jP(D  =  I  »  r< E X T  I 

t 

!  INvERT  MATRIX 
PGR  IR=1  TO  NP 

t 

!  FIND  T.dE  PIVOT  ELEMENT 

t 

P ;  VO  T  =  O 

FOR  I*IR  TO  NP  *  FDR  J=IR  10  NP  $  K*KP<I)  '»  L=Jp(J) 

I  tS  i  =Hbb  ( 3N l N  , L  )  ) -HbD  (r'iVu  i  > 

IF  TES i <  =  0  Then  GOTO  2220 
IF  TEST >0  THEN  GOTO  2190 
GOTO  2220 
IK*  I 

JK  =  J 

P IV0T*SN<  K , L  ) 

NEXT  J  4  NEXT  I 

♦ 

i  THE  PIVOT  ELEMENT  IS  SNU.L) 

'■  UPDATE  PIVOT  VECTORS 

K-KPUK)  ii  KP(IK)»KP(IR)  g  KPcIR)=K  L-JP  JK)  i  JP( jK > - JP ( IR >  »  JP(IS> 


■  '-UHPIJTE  NEW  ELEMENTS  uF  PIVOTAL  ROW 
,  F0R  TO  NCOL  9  SN<K,J)-SN<K, Jl/PIVUT  j)  1)0X1  J 
1  COMPOTE  REMAINING  ELEMENTS  OF  IR'TH  STEP 


2350  l 

236o  For  r-i  To  nP 

2370  IF  I  - K  =  0  THEN  GOTO  2440 

2380  AIL-SN( I ,L> 

2390  SN<I,L>  — (AIL/PIVOT) 

2400  FOR  J« I  TO  NCOL 

2410  IF  J-L-0  THEN  GOTO  2430 

2  A20  SN<I, J)-SNCI ,J>-AIL*SN(K, J> 

2430  NEXT  J 

2440  NEXT  I 

2450  ! 

2460  ! 

2470  !  INVERT  PIVOT  ELEMENT 

L48G  ! 

2490  $N(<  ,L)  =  I  /PIVOT  s>  NEXT  IR 

2500  RETURN 

2510  Hpiot : 

2520  !  . . 

2530  ’  • 

2540  ?  *  plotting  subroutine 

2550  !  . 

2560  ! 

2570  NuT= 1 

2580  CLnAR  ?  IF  NUMPlT-J  T~EN  GO . 0  ^7dU 
2530  ALPHA  6  DISP  "PARAMETERS  OF  FIT"  6  DISP 
2bOG  FOR  1=1  TO  NP  4  DISP  AA<I)  4  NEXT  I 
2b  1  0  IF  NSM  =  2  Then  GOTO  2b  4  l 
2b20  FOR  1=1  TO  NP  4  CUuL ( I ) =  AA<  I )  4  NEXT  I 
2d JO  GOTO  2650 

2b4u  FOR  1=1  TO  NP  4  HEAT ( I )  =AA(  I )  4  NEXT  I 
2650  CHI2=U 

2660  IF  NTEST<?  ThEN  GO'O  2700 
2670  SAR=4.2*<HEAT( 1 )-C3D^<  1  )  ) 

2660  IF  NS=2  THEN  SAR=2«SAR«h7 IME < JTURN ) +4 . 2- ( hEAT ( 2 > -COOL < 2 > > 

2b90  DISP  FNIVJDSC'SAR  IS"  )  ,  SAR  .  "Hat  ts  per  gran" 

2700  H OR  1=1  TO  N  4  Y2=0  4  PGR  j=1  IQ  NS+1 
2710  IF  NS  +  1 - J  =  0  AND  BT IME  C I >  =0  THEN  GOTO  2730 
2720  Y2  =  AA  C  J )  *BT  IfiE  ( I  >  (NS+1-D  +  Y2  4  GOTO  2740 
2/30  Y2=AA ( J )  +  Y2 
2/40  NEXT  J 

2750  CH 12= (67EMP ( I ) - Y2 )  2/Y2  +  CHI2 
2760  NEXT  I 

2770  DISP  "OBSERVED  CHI2=  ".CHI2,"  WITH  ".N+NS+I,"  D.F." 

2/dU  DISP  "DO  YUU  MANT  A  HARD  COPY  PLul'"  Y(tb>  Urt  N(U>' 

2790  INPUT  NYS4  IF  NUHSHO-1  THEN  GOTO  2890 

2600  DISP  "TIME  WILL  BE  PLOTTED  FROM  " , XMIN , " TO" , XMAX 

2610  DISP  "TEMPERATURE  WILL  BE  PLOTTED  FROM"  ,  YtllN  ,  "TO"  ,  YMAX 

2d20  DISP  "IF  YOU  ARE  SATISFIED  PRESS  KEiUKN" 

2830  INPUT  NYL IMS 

2840  IF  NYLIMS-""  THEN  GOTO  2890 

2850  DISP  "ENTER  TIME  RANGE  TO  BE  PLOTTED,  TNIN.TMAX” 

^bbU  1NHUI  XMIN, XMAX 

2870  DISP  "ENTER  TEMPERATURE  RANGE  TO  BE  PLOTTED" 

288U  INPUT  YMIN.YMAX 

2890  IF  NTEST>-  1  THEN  GOTO  3050 

2900  CLEAR 

2310  DISP  "CALCULATING  PLOT  SCALING  PARAMETERS  MaFNIVIDS< "PLEASE  WAIT") 
332U  AXMAX'lN;  (lGT  (  XMAX  >  ) 

2930  AYMAX-INT  ( LGT  (YMAX)) 

2340  Xr  *10' AXMAX  *  YF-10  AYMAX 

2350  XPWR4  =  VAL3  (AXMAX)  4  YPWRS'VAL's  (AYMAX) 

29no  *  -ii  ».  y+';n-l 


2970  IF  XMAX/XF>-  1  THEN  DELX-.2 
2980  IF  XMAX/XF>«  3  THEN  DELX-.5 
2990  IF  XMAX/XF>-  5  THEN  DELX-1 
3000  IF  YMAX/YF>-  1  THEN  DELY-.2 
3010  IF  YMAX/YF>-  3  THEN  DELT-.5 
3020  IF  YMAX/YF>-  5  THEN  DELY-1 
3030  XNUN-CEIL  ( CXMAX-XHIN)/(XF«DELX) > 

3040  YNUM-CEIL  ( < YMAX- YMIN) /< YF*DELY > > 

3050  IF  NYS-"N"  THEN  GOTO  3090 
3060  PLOTTER  IS  505 

3070  GRAPH  $  IF  NUHPLT-0  THEN  GCLEAR  @  DEG 
3U80  GOTO  3130 
3u90  PLOTTER  IS  1 
3100  GRAPH 

3110  !  IF  NUMPLT*0  AND  NTEST>0  THEN  GOTO  3110 

3120  IF  NUMPLT-0  THEN  GCLEAR 

3130  LOCATE  30, RATIO  *100-10,20,95 

3140  !  IF  NUMPlT*0  AND  NTEST>0  THEN  GOTO  3110 

3150  SCALE  XMIN/XF,XMlN/XF+XNUH*DELX, YMIN/YF, YMIN/YF  «-YNUM*DELY 

3lbU  !(•  N  I  LSI  >  0  ANU  NUMPlI^U  ImEN  bU  !  li  /3e>u 

3170  FXD  1 

3180  LGRID  DElX,DELY,XMIN/XF,YMIN/YF 
3190  AXS=(XMAX-XMlN)/(XF*iO> 

3/00  AYS=(YMAX-YMIN)/(YF*10) 

3210  MOVE  XMIN/XF+4.5*AXS,YMIN/YF-2*AYS  &  LORD  b 
3220  LABEL  USING  "K"  ;  "TIME  (SECS)  x  lE"iXPwRS 
J23U  HUVt  XMIN/X(--2*AXS,YMlN/YF+A.D*Atb 
3/AO  LORG  A  &  DEG  «  LDIR  90 

3250  LABEL  USING  "K"  ;  "TEMP . (DEGS  C>  x  1 E"A YPWR4 
3260  LDIR  0 

3270  IF  NPRQBE  =2  THEN  GOTO  33A0 
3280  LINE  TYPE  1 

3290  MOVE  ATInE< JMIN)/XF ,ATEMP< JM IN) 7 YF 
33UU  FOR  I-JMIN  TO  JMAX 
3/10  DRAW  ATIME<I)/XF,ATEMP(I>/YF 
3320  NEXT  I 

3330  IF  NPROBE-1  THEM  GOTO  3370 
33A0  LINE  TYPE  5 

3350  MOVE  BTIME(KMIN)/XF,BTEMP(KMIN)/YF 

3360  FOR  1-KM1N  TO  KMAX  *  DRAW  BTIME  ( I  )/XF  ,BTEMP(  I  >  /  YF  t  NEXT 
3370  IF  NYSLS«"N"  OR  NYSL$-"NQ"  THEN  GOTO  3A30 
3380  LINE  TYPE  A 

3390  MOVE  BTIME (801 >/XF ,BTEHP(80! )/YF 

3A00  FOR  1-801  TO  1000 

34 1 0  DRAW  BTIME( I) /XF ,BTEMP< I ) /YF 

3A20  NEXT  I 

3430  PEN  UP 

34A0  MOVE  XMIN/XF, YMIN/YF 

3450  IF  SAR-0  THEN  GOTO  3490 

3460  MOVE  XMIN/XF+4.5«AXS, YMAX/YF+. 1*AYS 

3470  LORG  4 

3480  LABEL  "SAR" ,SAR , "Watts  per  gran" 

3490  NTEST-NTEST+1 
3500  PAUSE 
3510  KEY  LABEL 
3520  RETURN 
3530  CHAIN  "DUMP" 

35au  END 
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INTRODUCTION  j 

Dielectric  heating  is  a  recognized  mechanism  for  the  induction  of  radiofrequency  (RF)  j 
bioeffects  considered  to  be  thermal  in  origin  •  that  is  attributable  to  temperature  increases. 
However,  apart  from  RF-heat  cell  killing  (e.g.  Sapareto  et  al.,1982,  Chang  et  al.,1987)  veryj 
little  information  is  available  on  quantitative  relationships  between  RF  bioeffects  at  the 
cellular  and  molecular  level  and  temperature  profiles  over  time  or  thermal  dosage  (TD).  The 
determination  of  RF  heating/cooling  curves  can  be  used  to  compute  the  specific  absorption 
rate  (SAR)  and  the  specific  absorption  (SA)  in  an  exposed  sample  (Stuchly  and  Stuchly, 
1986).  Biological  variables  and  exposure  conditions  can  be  controlled  in  experiments  with  in1 
vi fro  systems  to  a  degree  not  achievable  in  vivo.  An  attractive  model  to  study  RF  bioeffects 
is  the  transformation  of  lymphocytes  in  vitro,  provided  the  biological  variables  can  be  related 
to  dosimetric  quantities  that  characterize  exposure  (Czerski,1975;  Budd  and  Czerski,1985). 
To  accomplish  this,  we  designed  an  exposure  system  with  provision  for  real-time  temperature 
monitoring  with  RF-field  non-perturbing  temperature  probes.  The  exposure  system,  which 
will  be  described  in  detail  later,  has  multiple  sample  chambers.  To  allow  on-line  thermometry 
and  dosimetry,  one  of  these  chambers  is  used  as  a  site  for  a  non-perturbing  temperature  probe. 
From  the  temperature  (T)/time(t)  history  of  the  sample  chamber,  an  exposure  dosage  can  be 
determined  and  described  in  terms  of  TD,  SAR,  S  A,  or  electric  field  strength.  The  system  can 
be  used  for  studies  of  RF-bioeffects  in  any  tissue  culture  cell  line  or  other  in  vitro  biological 
sample.  The  monitoring  component  can  be  applied  to  the  study  of  temperature-dependent 
effects,  irrespective  of  the  modality  used  for  heating. 

EXPOSURE  SYSTEM 

A  schematic  diagram  of  the  system  for  microwave  exposure  of  samples  is  shown  in  Fig.  1. 


The  output  from  *  CW  Hewlett-Packard  Model  8616A  oscillator  feeds  a  Hughes  Model  1177H 
TWT  amplifier.  Pulsed  or  amplitude-modulated  waveforms  are  obtained  through  the  use  of 
a  Hewlett-Packard  Model  8403 A  p-i-n  modulator  and  driver.  For  high  peak  power,  but  low 
duty-cycle,  pulsed  exposures  the  signal  can  be  further  amplified  with  a  lkW  amplifier  (MCL, 
La  Grange,  Ill.)  The  amplified  signal  enters  a  aborted  section  of  S-band  rectangular  waveguide 
300  mm  long  through  a  matched  coaxial  feed.  The  waveguide  sample  holder  receives  its 
power  through  an  isolator  (not  shown  specifically  in  Fig.  1)  and  a  dual-directional  coupler 
that  allows  forward  and  reflected  powers  to  be  measured  with  Hewlett-Packard  Model  432A 
power  meters,  and  waveforms  to  be  monitored.  A  sample  holder  for  biological  specimens 
is  supported  centrally  in  the  waveguide  in  a  block  of  low-density  polystyrene  that  fills  the 
cross-section  of  the  guide.  An  identical  waveguide  section  with  an  identical  plastic  sample 
holder  is  used  for  sham  exposures.  During  microwave  exposures,  both  waveguide  assemblies 
were  housed  in  a  C02  tissue  culture  incubator  thermostated  at  S7.0°C.  Various  sample 
holders  were  tested:  a  4-chamber  plastic  tissue  culture  slide  (Miles  Scientific  #  4804)  was 
found  satisfactory  in  our  application.  Each  chamber  holds  1  ml.  A  miniature  thermistor 
temperature  probe  enters  each  waveguide  section  through  a  small  hole  in  the  top  shorting 
plate.  The  hole  is  drilled  near  the  guide  wall,  at  the  center  of  the  shorter  dimension  -  this  is 
a  low  electric  field  point  and  the  hole  produces  minimal  disturbance.  The  active  end  of  the 
thermistor  probe  enters  the  sample  in  one  of  the  4  chambers  of  the  tissue  culture  dish.  We 
could  move  the  probe  from  chamber  to  chamber  to  check  exposure  uniformity.  Two  types  of 
non-field  perturbing  probe  have  been  used  in  this  way:  A  Narda  Model  801  IB  non-perturbing 
double  temperature  probe  or  two  Vitek  Model  101  probes.  The  temperatures  in  both  the 
exposed  and  sham-exposed  samples  are  recorded  continuously  during  an  experimental  run. 
The  temperature  probes  are  connected  to  a  Hewlett-Packard  relay  activator,  Model  59306A, 
which  is  itself  connected  to  a  Keithley  Model  192  DVM.  Both  the  relay  activator  and  DVM 
are  under  the  control  of  a  Hewlett-Packard  Model  86  desk-top  computer  through  the  HPIB 
(IEEE-488  bus)(Fig.2). 


Fig. 2.  Schematic  diagram  of  multiple-chamber  sample  bolder  with  son- 
perturbing  temperature  measurement  and  automated  dosimetry  system. 


Under  computer  control,  the  temperatures  of  exposed  and  sham-exposed  samples  are 
recorded  sequentially  and  stored  in  memory.  Since  the  temperature  of  sham-exposed  samples 
is  essentially  constant,  it  is  monitored  less  frequently  than  the  temperature  of  the  exposed 
samples.  A  typical  experimental  protocol  would  involve  10  temperature  readings  of  the  ex¬ 
posed  sample  for  every  one  reading  of  the  sham-exposed.  For  conventional  thermal  exposures, 
the  temperature  can  be  elevated  by  increasing  the  temperature  of  the  incubator. 

TEMPERATURE/TIME  ANALYSIS 

During  a  typical  experiment,  the  temperature  is  recorded  at  regular  intervals,  with  a 
minimum  measurement  interval  <  Is,  before  the  beginning  of,  during,  and  after  the  exposure. 
The  “on"  and  “off"  times  of  exposure  are  recorded  on  the  computer  by  the  operator  using  a 
“soft-key"  interrupt  capability.  At  the  conclusion  of  a  run,  the  T(t)  behavior  is  analysed  to 
determine  SAR 

If  the  SAR  is  S  (Wkg-1)  and  the  specific  heat  of  the  sample  is  C(Jkg“IK_1),  then  the 
SAR  can  be  determined  from  the  change  in  heating  rate  when  microwave  power  is  applied. 


In  practice,  this  analysis  is  performed  automatically.  At  the  end  of  a  run,  the  point  of 
inflection,  or  turning  point,  in  T(t)  is  found  numerically.  The  temperature  profile  to  the  left 
and  right  of  this  point  is  fitted  by  least-squares  to  a  linear  or  quadratic  function  and  the 
change  in  slope  at  the  beginning  of  the  exposure  period  yields  the  SAR.  Figures  3  and  4  are 
examples  of  such  a  procedure. 
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Fig. 3.  Temperature/time  history  of  an 
exposed  sample  showing  linear  fits  to  cool¬ 
ing  and  heating  portions  of  curve  for  dosime 
try. 


Fig. 4.  Temperature/time  history  of  ex¬ 
posed  (upper  curve)  and  sham-exposed 
samples  (lower  curve).  The  dosimetry 
for  the  exposed  sample  has  been  deter¬ 
mined  from  a  quadratic  fit  to  the  heating 
curve  during  microwave  exposure  and  a 
linear  fit  to  the  equilibrated  portion  of 
the  curve  prior  to  exposure. 


Figure  3  shows  the  temperature/time  history  of  a  sample  that  was  cooling  prior  to  the  start 
of  microwave  exposure.  The  linear  fits  to  the  cooling  and  heating  portions  of  the  curve  yield 
the  SAR.  Figure  4.  shows  the  heating  of  a  sampled  that  was  quite  well  equilibrated  befcre 
exposure.  Because  the  sigual/noise  ratio  is  high  the  SAR  was  determined  by  a  quadratic 
fit  to  the  heating  portion  of  the  curve.  The  quality  of  the  fit  sear  the  turning  point  can 
always  be  examined  to  determine  whether  the  slope  value  is  realistic.  However, unless  the  , 
temperature/time  profile  has  very  low  noise,  it  is  generally  better  to  use  the  linear  fit  to  T|t) 
near  the  turning  point  or  a  biased  estimate  of  the  slope  can  result.  Figure  4  also  shows  the 
temperature/time  history  of  the  sham-exposed  sample  to  illustrate  its  temperature  stability. 

RESULTS  AND  DISCUSSION 

The  advantages  of  the  procedure  described  above  are  severalfold.  The  sample  need  not 
be  equilibrated  before  SAR  is  determined.  Repeated  determination  of  SAR  can  easily  aid 
quickly  be  made,  which  allows  the  SAR  uniformity  from  one  sample  chamber  to  another 
to  be  determined.  Exact  knowledge  of  the  microwave  power  is  not  required:  reproducible 
exposures  can  be  made  at  known  SAR  by  using  the  measured  forward  power  corresponding 
to  a  given  SAR  measurement,  provided  the  experimental  arrangement  is  not  altered  between 
exposures.  Dosimetry  is  not  affected  by  other  losses  in  the  system.  The  waveguide  exposure  ' 
system  with  four  sample  chambers  was  found  to  give  an  SAR  uniformity  from  one  chamber 
to  the  other  within  10%.  Thus,  experimental  samples  can  be  exposed  and  examined  in 
triplicate,  the  fourth  chamber  being  used  for  the  insertion  of  the  temperature  probe. 

The  system  is  biocompatible  and,  depending  on  cell  line,  cell  density,  and  medium,  allows 
continuous  exposures  of  several  days  duration.  The  biocompatibility  of  this  arrangement  was 
proven  by  studying  the  growth  of  human  lymphocytes  under  various  conditions  within  the 
waveguide.  As  well  as  providing  biocompatible  exposure  conditions,  the  system  is  flexible: 
it  allows  the  exposure  of  tissue  culture  cells  growing  in  suspension  or  in  monolayers,  and  can 
easily  be  used  with  cell-free  samples.  The  exposure  system  is  relatively  simple,  and  can  be 
assembled  from  off-the-shelf  components.  Its  capabilities  can  be  expanded  by  introducing 
computer  control  of  RF  power  input  b<\sed  on  feedback  from  temperature  measurements. 
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ABSTRACT 


There  ate  contradictory  reports  In  the  literature  on  the  effects  of  mlcro- 
wa»e  exposure  on  transformation  of  lymphocytes.  To  resolve  this  question, 
human  lymphocytes  have  been  exposed  1  n  vitro  to  2450  WHx  radiation  under 
carefully  controlled  conditions  at  various  power  levels  and  for  periods  up 
to  120  hour  s.  Lymphocyte  tr  information  In  exposed  samples  was  compared  to 
control  samples,  sham-exposed  or  Incubated  In  a  conventional  tissue  culture 
Incubator.  Heparinized  peripheral  blood  samples  were  obtained  by  venipunc¬ 
ture  from  healthy  volunteers,  lymphocytes  were  Separated  by  gredlent  cen¬ 
trifugation  (Lymphocyte  Separation  Medium,  Litton  Blonetlcs)  and  Incubated 
In  chromosome  medium  1A  without  phytohemagglutinin  (Glbeo)  et  a  Standard 
concentration  of  10  cells/ml  medium.  Using  the  automated  dosimetry  micro- 
wave  exposure  system  described  In  e  companion  paper  (Joyner  et  el.,  this 
Congress)  the  temperature  profile  over  tine  was  recorded  for  each  sample, 
and  average  end  peak  specific  absorption  rates  (SAft),  specific  absorption 
(SA),  end  thermal  dosage  (TO)  were  calculated.  At  harvest,  cells  were 
counted  and  their  viability  was  tested  using  neutral  red,  Janus  green,  and 
trypan  blue  staining.  Standardized  cytologlcal  preparation  were  made  using 
a  "Cytospln  11"  (Shandon)  cytocentr Huge,  alt-dried,  fixed  in  methanol  and 
stained  with  a  combination  of  Wright  and  Glemsa  stains.  The  percentages  of 
untransformed  (small)  lymphocytes.  Intermediate  and  lymphoblastold  cells 
wre  determined,  based  on  morphological  features,  Quantitated  using  an 
*Opto»ex*  Image  analyzer.  The  correlation  of  lymphocyte  transformation 
with  SAft ,  sa,  and  TO  was  examined.  The  Implications  of  the  results  for  a 
‘thermal*  vt  "nonthermal*  mechanism  of  interaction  will  be  discussed. 
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The  effects  of  microwaves  on  the  acetylcholine  (ACh)  receptor  channel  have 
been  examined  recording  ACh-induced  currents.  Experiments  were  performed  on 
single  3-5  day  old  unstriated  and  non -innervated  quail  myotubes  .  The 
properties  of  single  ACh  activated  channels  were  investigated  using  the 
patch-clamp  technique  with  high  resistance  seals  between  the  membrane  and 
the  pipette  in  cell  attached  configuration.  Sylgard  coated  pipette.*  were 
filled  with  acetylcholine  at  the  concentration  5x10  aM.  Single  channel  w#r« 
analyzed  by  the  threshold  crossing  method  Of  microwaves  (7  itH  it  10.75  GHz ) 
were  applied  using  a  horn  antenna  (16  dB  gain)  at  a  distance  of  25  cm.  from 
the  sample.  30*  inclinated  to  the  ground  and  with  the  E  vector  lying  on  the 
plane  of  incidence  In  the  control,  before  the  exposure,  ACh-activated 
channel  current  of  1  74  i  0.01  pA  (mean  i:  a.e.m.;  740  events).  a  mean 
channel  open  time  of  3.46  i  0.13  ms,  a  mean  closad  time  of  18.7  i  1.3  ms 
were  measured  When  myotubes  were  exposed  to  toi  for  120  s,  the  closed  time 
increased  to  50 . 1  1  5.1  ms  (300)  within  90s  from  the  onset.  Channel  open 
time  (3.43  +  0.16  ms)  and  current  (1.77  it  0.01  pA)  were  not  affected  by 
irradiation  After  the  end  of  exposure  a  partial  recovery  of  the  activity 
was  obtained  It  is  concluded  that  exposure  to  mi c reserves  increases  the  rate 
of  desensitization  and  decreases  the  channel  opening  probability.  Other 
experiments  indicate  that  this  effect  develops  only  within  range  of  field 
intensi ty . 
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B-CP-f  EFFECTS  OF  2.45  GHz  EXPOSURE  ON  SPONTANEOUS  HUMAN  LYMPHOCYTE 

TRANSFORMATION  IN  VITRO.  E  C  Elson,  Walter  Reed  AIR,  Silver  Spring.  MD 

20905  E  M  Czerske,  C  C.  Davie,  Electrical  Engineering  Dept.,  University 
of  Maryland.,  College  Perk.  K3  20742,  P.  Ctereki.  CDRH/FDA.  Rockville,  MD 
20857. 

Him  an  peripheral  blood  lymphocytes  in  chromoeome  medlun  1A  without  R1A 
(Gibco)  have  been  exposed  to  continuous  (CW)  or  pulsed  (BJ)  2.  *5  GHr 

microwaves  using  tha  automated  dosimetry  exposure  system  described  in  a 
ccxnpanion  paper  (Joyner  et  al . ,  this  meeting)  at  heating  mi  nen-heating  SAR 
levels  for  periods  up  to  120  hours.  At  harvest,  standardized  cytologlcel 
preparations  were  mede  using  a  cytospin  II  cytocentrifuge,  air-dried,  fixed 
in  methanol  and  stained  with  Wright/Giemsa  stain.  The  percentages  of 
untrana formed  lymphocytes,  intermediate  and  transformed  lynyhob las told  cells 
were  determined  based  on  morphological  features  qusntltated  using  an 

“Optomax"  image  analyser  and  coopered  b  a  tween  exposed  (CW  or  fW)  and  ah  art - 
exposed  samples  W  or  IV  exposure  at  the  ease  average  SAR  resulting  In 
•  imilar  temperature  profiles  over  time  affect*  lymphocyte  transformation 
differently  Mechanistic  implications  of  this  finding  are  discussed. 
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PW-C5  AUTOMATED  DOSIMETRY  OF  MICROWAVE  AND  THERMAL  BCPOSURE  OF -BIOLOGI¬ 
CAL  SAMPLES  IN  VITRO.  K.H.  Joyr.tr,  TtUcom,  Melbourne,  Australia.  E.C. 
Elton,  Walter  Rt«d  AIR,  Silvtr  Spring,  H>  20 905.  C.C.  Dmr It,  E.M.  Cxerike,* 
Electrical  Eng  Dtpt.,  U  of  MD.,  College  Park,  M)  20742,  P.  Cserski. 
CDRH/FDA,  RockviUt.  tf)  20657, 

A  microwave  exposure  and  dosimatric  system  will  bt  daacrlbad  that  provides 
non-perturbing  automated  recording  of  the  teqperature  history  of  exposed  «rd 
control  samples  in  vitro.  Records  of  the  temperature (T) /time (t )  behavior  of 
both  sample  and  control  are  made  before,  during,  and  after  exposure.  From 
the  recorded  T(t)  versus  t  during  heating  and  cooling  the  specific  absorp¬ 
tion  rate  (SAR)  for  microwaves  (or  the  equivalent  energy  abaorption  during 
conventional  heating)  can  be  determined  The  exposure  systen  oaisista  of  two 
shorted  sections  of  S-bend  rectangular  waveguide  with  matched  coaxial  feeds 
placed  In  a  C02  tissue  culture  incubator  thermostated  at  37.5*C.  Sample 
holders  are  supported  centrally  in  the  waveguide  in  blocks  of  low-density 
expanded  polystyrene.  An  oscillator,  TWT  amplifier,  Isolator,  and  forward 
and  reflected  power  meters,  provide  the  power  to  the  exposure  waveguide.  The 
second  waveguide  houses  the  control  samples.  Temperatures  in  one  sample 
chamber  In  each  waveguide  ere  recorded  with  e  Narda  ttxiel  801  LB  ncn -perturb¬ 
ing  probe  The  temperature  readings  ere  digitised  and  stored  in  a  desktop 
computer.  Subsequent  linear  least  squares  v^lysis  of  T(t)  versus  t  with  the 
assumption  of  constant  power  absorption  aid  Newton's  La#  of  cooling  provides 
a  very  reliable  value  for  the  SAR  On-line  determination  of  the  SAR  in  this 
way  allows  accurate  assessment  of  exposure  uniformity  from  one  sample 
chamber  to  the  other,  end  permits  accurate  calibration  of  the  system  in 
terms  of  e  SAR.' forward-power  ratio. 
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PW-Ce  INTENSITY  DEPENDENCE  OF  THE  INDUCTION  OF  ACETYLCHOLINE  ESTERASE 
ACTIVITY  BY  MODULATED  RADIOFREQUENCY  RADIATION  S.  K.  Dutte,  B.  Ghosh* ,  C 
F  Blackman  .  Botany  Department  and  Cancer  Research  Center.  Howard  Univer¬ 
sity.  and  ^Health  Effects  Research  Laboratory.  USEPA,  Res.  Tn.  Pk .  ,  NC . 

Radio  frequency  radiation  (RFR)  at  915  and  at  147  MHz,  when  sinusoidally 
amplitude  modulated  (AM)  at  16  Hz,  has  been  shown  to  nhaiced  the  release  of 
calcium  ions  front  neurobias  tome  cells  in  culture.  The  dose  response  for  this 
effect  is  highly  unusual,  consisting  of  two  power  droity  ‘Windows”  in  ehich 
enhanced  efflux  occurs,  separated  by  power  density  regions  where  no  effect 
is  observed  A  similar  result  has  also  bean  reported  for  s  chicken  brain- 
tissue  preparation  In  order  to  explore  ths  physiological  importance  of 
these  findings,  we  have  examined  the  impact  of  exposure  on  a  membrane- bound 
enzyme,  acetylcholine  esterase  ( AChl ) .  which  is  intimately  involved  with  the 
acetylcholine  neurotransmi tter  system  Neuroblastoma  cells  (K3-106),  exposed 
for  30-minutes  to  14  7-ffjz  radiation,  AM  at  16  Hz,  damnstrated  mtunced  AChE 
activity,  as  assayed  by  a  procedure  using  ^C* labeled  acetylcholine.  This 
enhanced  activity  was  observed  only  within  a  narrow  time  window  aaverel 
hours  after  the  cells  were  Beaded ,  and  only  when  the  exposure  occurred  at 
power  densities  that  ware  affective  for  calcium  ions.  This  result  indicates 
that  RFR  can  perturb  nervous-systero-der ived  calls  in  culture  to  affect  both 
calcium  ion  release  and  AChE  activity  in  a  cocrooo  doaa- dependent  manner. 
Both  the  physiological  significance  and  the  underlying  mechanism  of  action 
of  these  low  intensity  fields  will  need  to  be  reviewed  in  light  of  these 
findings  We  acknowledge  support  from  an  Institutional  Training  Grant  from 
Howard  University,  Cooperative  Agreement  CR812100  from  IEEPA,  and  Inter eg «n- 
cy  Agreement  DE-A106-67RL11374  from  DOE. 
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AN  AUTOMATED  DOSIMETRY  SYSTEM  FOR  MICROWAVE  AND 
THERMAL  EXPOSURE  OF  BIOLOGICAL  SAMPLES  IN  VITRO 
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Abstract— A  waveguide  exposure  system  with  automated  sample  temperature  measoreuieot  la  described  This 
system  provides  on-line  determination  of  tbe  temperature  profile  over  time  of  biological  samples  fa  ritro.  It  allows 
automated  computation  of  the  specific  absorption  rate  determined  from  beating  /cooling  curves,  uses  miaimally- 
perturbing  thermometry,  b  biocompatiblr  and  can  be  used  for  measurements  of  both  microwave  aad  conventional 
beating 


INTRODUCTION 

Dielectric  heating  is  a  recognized  mechanism  for  the 
induction  of  radiofrequency  (RF)  bioeffects  considered 
to  be  thermal  in  ongin — that  is,  attributable  to  temper¬ 
ature  increases  However,  apart  from  RF-heat  cell  killing 
(Saparcto  et  a]  1982,  Chang  el  al.  1987),  very  little  in¬ 
formation  is  available  on  quantitative  relationships  be¬ 
tween  RF  bioeffecis  at  the  cellular  and  molecular  level 
and  the  temperature  profile  over  time  (TPT).  The  deter¬ 
mination  of  RF  healing /cooling  curves  can  be  used  to 
compute  the  specific  absorption  rate  (SAR  W  kg'1)  and 
the  specific  absorption  (SA  J  kg'1)  in  an  exposed  sample 
(Stuchly  and  Stuchly  1986).  Biological  variables  and  ex¬ 
posure  conditions  can  be  controlled  in  experiments  with 
in  vitro  systems  to  a  degree  not  achievable  in  vivo.  An 
attractive  model  to  study  RF  bioeffects  is  the  transfor¬ 
mation  of  lymphocytes  in  vitro,  provided  the  biological 
variables  can  be  related  to  dosimetric  quantities  that 
characterize  exposure  (Czerski  1975;  Budd  and  Czerski 
1985 ).  To  accomplish  reliable  dosimetry,  we  designed  an 
exposure  system  with  provision  for  real-time  temperature 
monitoring  with  RF-field  non-perturbing  temperature 
probes.  The  exposure  system,  which  will  be  described  in 
detail  later,  has  multiple  sample  chambers.  To  allow  on¬ 
line  thermometry  and  dosimetry,  one  of  these  chambers 
It  used  as  a  site  for  a  non-perturbing  temperature  probe. 


From  the  temperature  ( T)/time(f )  history  of  the  sample 
chamber,  an  exposure  dosage  can  be  determined  and  de- 
acribed  in  terms  of  TPT,  SAR,  SA,  or  electric  field 
strength.  Tbe  system  can  be  used  for  studies  of  RF  bioef- 
fects  in  any  tissue  culture  cell  line  or  other  in  vitro  bio¬ 
logical  samples.  The  temperature  monitoring  pan  of  the 
system  and  the  associated  automated  data  analysis  can 
be  applied  to  the  study  of  temperature-dependent  effects, 
irrespective  of  the  modality  used  for  heating. 

EXPOSURE  SYSTEM 

A  schematic  diagram  of  the  system  for  microwave 
exposure  of  samples  is  shown  in  Fig.  1 .  Tbe  output  from 
a  CW  oscillator'  feeds  a  tnvelling-wave-tube  ("FAT )  am¬ 
plified.  Pulsed  or  amplitude-modulated  waveforms  are 
obtained  through  the  use  of  a  p-i-n  modulator  and  driver# 
For  high  peak  power,  but  low  duty-cycle,  pulsed  exposures 
the  signal  can  be  further  amplified  with  a  I  kW  ampli¬ 
fier*  *.  The  amplified  signal  enters  a  shoned  section  of  S- 
band  rectangular  waveguide  300  mm  long  through  a 
matched  coaxial  feed.  The  waveguide  sample  holder  re¬ 
ceives  its  power  through  an  isolator  and  a  dual -directional 
coupler  that  allows  forward  and  reflected  powers  to  be 
measured  with  power  meterstt  aad  waveforms  to  be 
monitored  with  a  diode  detector#  and  oscilloscope  A 
ample  bolder  for  biological  yeomens  is  supported  cew 
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Fig  1  Schematic  diagram  of  essentia]  components  of  microwave 
exposure  system 


Microwave  Powe'  In 


Fig  2  Waveguide  exposure  assembly  showing  location  of  four- 
chamber  sample  holder 


trally  in  the  waveguide  in  a  block  of  low-density  polysty¬ 
rene  that  fills  the  cross-section  of  the  guide,  as  shown  in 
Fig.  2.  An  identical  waveguide  section  with  an  identical 
plastic  ample  holder  is  used  for  (ham  exposures  During 
microwave  exposures,  both  waveguide  assemblies  were 
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Fig  3  Schematic  diagram  of  four-chamber  sample  holder  with 
minimally-perturbing  temperature  measurement  and  automated 
dosimetry  system 


boused  in  a  CO;  tissue  culture  incubator  thermostated  at 
37.0°C.  Various  ample  holders  were  tested  a  four- 
chamber  plastic  tissue  culture  slide§§.  shown  in  detail  in 
Fig  3.  was  found  atisfactory  in  our  application.  Each 
chamber  holds  I  mL.  Since  this  ample  holder  occupies 
only  about  10*  of  the  cross-sectional  area  of  the  wave¬ 
guide,  field  uniformity  from  one  chamber  of  the  ample 
holder  to  the  other  ihould  be  within  about  ±6*  of  the 
average  and  power  densities  should  be  within  ±12*  of 
the  average  This  expected  exposure  uniformity  was  borne 
out  by  quantitative  dosimetry  A  miniature  thermistor 
temperature  probe  enters  each  waveguide  section  through 
a  small  hole  in  the  top  shorting  plate.  The  hole  is  drilled 
near  the  guide  wall  at  the  center  of  the  shorter  dimension 
This  is  a  low  electric  field  point,  and  the  hole  produces 
minimal  disturbance.  The  active  end  of  the  thermistor 
probe  enters  the  ample  in  one  of  the  four  chambers  of 
the  tissue  culture  dish.  We  could  move  the  probe  from 
chamber  to  chamber  to  check  exposure  uniformity  Two 
types  of  minimally-field-perturbing  probe  have  been  used 
in  this  way:  a  Narda  Model  801  IB  non-perturbing  double 
temperature  probe’  or  two  Vitek  Model  101  probesTT. 
The  temperatures  in  both  the  exposed  and  sham-exposed 
samples  are  recorded  continuously  during  an  experimental 
run  The  temperature  probes  are  connected  to  a  relay  ac¬ 
tivator#**  which  is  connected  to  a  digital  voltmeter 
(DVM  )•  •  *.  Both  the  relay  activator  and  DVM  are  under 
the  control  of  •  desk -top  computerf+t  through  the  IEEE- 
488  bus  Under  computer  control,  the  temperatures  of 
exposed  and  sham-exposed  ramples  are  recorded  sequen¬ 
tially  and  stored  in  memory.  Since  the  temperature  of 
iham-exposed  umples  is  essentially  constant,  it  is  mon- 
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itored  less  frequently  than  the  temperature  of  the  exposed 
samples  A  typical  experimental  protocol  would  involve 
10  temperature  readings  of  the  exposed  sample  for  every 
one  reading  of  the  sham-exposed  For  conventional  ther¬ 
mal  exposures,  the  temperature  can  be  elevated  by  in¬ 
creasing  the  temperature  of  the  incubator. 

TEMPERATURE /TIME  ANALYSIS 

During  a  typical  experiment,  the  temperature  is  re¬ 
corded  at  regular  intervals,  with  a  minimum  measurement 
interval  of  <1  s  before,  during  and  after  the  exposure. 
The  ‘‘on’*  and  "off"  times  of  exposure  art  recorded  on 
the  computer  by  the  operator  using  a  “soft-key"  interrupt 
capability.  At  the  conclusion  of  a  run,  the  T(t)  behavior 
is  analyzed  to  determine  SAR 

If  the  SAR  is  S(W  kg' ')and  the  specific  heat  of  the 
sample  is  C  ( J  kg  * '  K. ' 1 ) ,  then  the  rate  of  heating  during 
microwave  exposure,  for  small  temperature  variations 
over  which  C  can  be  assumed  to  be  constant,  is 


If,  in  the  absence  of  deliberate  heaung.  a  sample  is  above 
or  below  its  equilibrium  temperature.  T with  its  sur¬ 


roundings,  then  for  small  temperature  differences  from 
equilibrium,  its  natural  rate  of  temperature  change,  from 
Newton's  law  of  cooling,  is  of  the  form 

where  a  is  a  constant  that  depends  on  sample  geometry, 
insulation,  etc.  In  general,  when  microwave  power  is  ap¬ 
plied  to  a  sample  that  was  not  in  thermal  equilibrium 
with  its  surroundings  at  the  start  of  exposure: 

dT  S 

(3) 

The  solution  to  eqn  (3)  can  be  written  in  the  form 

T  -  T0  «  (~  +  T*,  -  7o)(l  -  e"*'),  (4) 

where  T0  is  the  control  temperature  at  t  «  0.  Either  eqn 
(3)  or  eqn  (4)  can  be  used  to  determine  the  SAR  From 
eqn  (3)  we  have 
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Time  (b)  x  100 

Fig  4.  Temperature /time  profile  of  a  sample  that  was  subjected  to  RF  heating  (solid  curve).  The  dashed  curves 
are  least-squares  linear  fits  to  the  experimental  data  To  the  left  of  the  time  at  which  RF  poweT  was  applied,  the 
sample  was  cooling  The  change  in  slope  at  the  initiation  of  the  RF  beating  provides  the  quantitative  dosimetry. 
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where  subscript  5  indicates  the  application  of  microwave 
power  corresponding  to  SAR  S.  Therefore,  determination 
of  the  rate  of  temperature  change  before  and  during  ap¬ 
plication  of  microwave  power,  or  during  and  after  the 
application  of  microwave  power,  determines  the  SAR. 

If  the  sample  is  in  thermal  equilibrium  before  ex¬ 
posure  begins,  eqn  (4)  reduces  to  the  simpler  form 

T- -e-').  (6) 

which  can  be  used  to  find  5  from  T(t).  In  general,  the 
approach  using  eqn  ( 5 )  is  preferable,  since  for  times  near 
I  ■  0,  the  exponential  can  be  expanded  in  quadratic  form. 

In  practice,  this  analysis  is  performed  automatically. 
At  the  end  of  a  run.  the  point  of  inflection,  or  turning 
point,  in  T(t)  is  found  numerically.  This  is  done  by 
checking  the  change  in  slope  of  the  temperature  as  a  func¬ 
tion  of  time  In  practice,  this  change  in  slope  occurs  im¬ 
mediately  upon  the  application  or  removal  of  microwave 
power,  within  the  time  resolution  of  the  data  acquisition 
system  The  temperature  profile  to  the  left  and  right  of 
this  point  is  fitted  by  least-squares  to  a  linear  or  quadratic 
function,  and  the  change  in  slope  at  the  beginning  of  the 
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exposure  period  yields  the  SAR.  Figures  4  and  5  are  ex¬ 
amples  of  such  a  procedure.  Figure  4  shows  the  temper¬ 
ature/time  history  of  a  sample  that  was  cooling  prior  to 
the  start  of  microwave  exposure.  The  linear  fits  to  the 
cooling  and  beating  portions  of  the  curve  yield  the  SAR. 
Figure  S  shows  the  heating  of  a  sample  that  was  quite  well 
equilibrated  before  exposure.  Because  the  signal  /noise 
ratio  is  high,  the  SAR  was  determined  by  a  quadratic  fit 
to  the  heating  portion  of  the  curve.  The  quality  of  the  fit 
near  the  turning  point  can  always  be  examined  to  deter¬ 
mine  whether  the  slope  value  is  realistic.  However,  unless 
the  temperature /time  profile  has  very  low  noise,  it  is  gen¬ 
erally  better  to  use  the  linear  fit  to  T(t)  near  the  turning 
point,  otherwise  a  biased  estimate  of  the  slope  can  result. 
Figure  5  also  shows  the  temperature /time  history  of  the 
sham-exposed  sample  to  illustrate  its  temperature  stability. 

RESULTS  AND  DISCUSSION 

The  advantages  of  the  procedure  described  above  are 
severalfold  The  sample  need  not  be  equilibrated  before 
SAR  is  determined  Repeated  determination  of  SAR  can 
easily  and  quickly  be  made,  which  allows  the  SAR  uni¬ 
formity  from  one  sample  chamber  to  another  be  deter- 
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Fig  5  Temperature /ume  profile  of  exposed  (upper  soIkJ  curve)  end  tham-expoeed  samples  (lower  dashed 
curve )  The  dosimetry  for  the  exposed  sample  has  been  determined  from  ■  quadratic  fit  to  the  beating  curve 
during  microwave  exposure  sad  a  linear  fit  to  the  equilibrated  portion  of  the  curve  prior  to  exposure  (dashed 

upper  curve ) 
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mined  Exact  knowledge  of  the  microwave  power  is  not 
required;  reproducible  exposures  can  be  made  at  known 
SAR  by  using  the  measured  forward  power  corresponding 
to  a  given  SAR  measurement,  provided  the  experimental 
arrangement  is  not  altered  between  exposures.  Dosimetry 
is  not  affected  by  other  losses  in  the  system.  The  waveguide 
exposure  system  with  four  sample  chambers  was  found 
to  give  an  SAR  uniformity  from  one  chamber  to  the  other 
within  10%.  Thus,  experimental  samples  can  be  exposed 
and  examined  in  triplicate,  the  fourth  chamber  being  used 
for  the  insertion  of  the  temperature  probe. 

The  system  is  biocompatible  and,  depending  on  cell 
line,  cell  density  and  medium,  allows  continuous  expo¬ 
sures  of  several  days  duration  The  biocompatibility  of 
this  arrangement  was  proven  by  studying  the  growth  of 
human  lymphocytes  under  various  conditions  within  the 
waveguide  A  concentration  of  10*  cell  mL'1  in  chro¬ 
mosome  medium  lAftt  with  or  without  phytohemag- 
glutin  was  incubated  for  72  or  120  h  in  the  exposure  sys¬ 
tem  with  no  applied  RF  power.  Cell  morphology,  lym- 
phoblastoid  transformation,  mitotic  index  and  cell 
viability  tested  by  neutral  red,  Janus  green  and  trypan 
blue  stains  did  not  differ  from  those  in  lymphocyte  cul¬ 
tures  from  the  same  donor  cultured  in  a  conventional 
CO;  incubator  according  to  the  protocol  provided  by 
Gibco  with  their  media  As  well  as  providing  biocompat¬ 
ible  exposure  conditions,  the  system  is  flexible.  It  allows 
the  exposure  of  tissue  culture  cells  growing  in  suspension 
or  in  monolayers  and  can  easily  be  used  with  cell-free 
samples 

The  exposure  system  is  relatively  simple  and  can  be 
assembled  from  off-the-shelf  components  Its  capabilities 
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can  be  expanded  by  introducing  computer  control  of  RF 
power  input  based  on  feedback  from  temperature  mea- 
gurements  in  a  similar  way  to  the  system  described  by 
Chang  et  al  ( 1987 ).  This  system  used  a  stripline  operating 
at  915  MHz  and  allowed  uniform  exposure  of  10  mL 
tissue  culture  samples.  The  temperature /time  profile  of 
beating  could  be  pre-programmed,  and  biocompatibility 
allowed  exposures  of  several  days  duration.  A.  W.  Guy 
(1977)  designed  a  system  for  exposure  of  5  mL  samples 
in  the  range  from  dc  to  1  GHz.  Field  strength  dosimetry 
was  obtained  from  on-line  measurement  of  feed-line 
impedance.  The  temperature  of  the  sample  could  be 
maintained  at  a  preset  level  using  a  circulating  liquid  heat 
exchanger.  No  data  on  biocompatibility  of  this  system 
were  provided,  however  the  information  given  suggests 
suitability  only  for  short-term  exposures.  J.  C.  Lin  ( !  9Tri  > 
and  Chen  and  Lin  (1978)  proposed  the  use  of  micropi¬ 
pettes  immersed  in  a  fluid-filled  waveguide  irradiation 
chamber  for  in  vitro  studies  under  controlled  temperature 
and  dosimetric  conditions.  This  system  was  also  only 
suitable  for  short-term  exposures.  J.  W.  Allis  et  al.  (1975 ) 
installed  a  waveguide  exposure  system  in  a  Cary  1 5  spec¬ 
trophotometer  allowing  simultaneous  exposure  and  spec¬ 
trophotometnc  biochemical  determinations.  Our  system, 
as  well  as  those  described  in  the  papers  quoted  here,  pro¬ 
vides  a  wide  range  of  possibilities  for  in  vitro  studies  under 
controlled  conditions,  including  precise  RF  dosimetry 
By  using  one  of  these  systems,  or  combining  desired  fea¬ 
tures  of  some  of  them,  most  requirements  for  the  study 
of  RF  bioeffects  at  the  cellular  and/or  molecular  level  can 
be  met. 
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